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There are many interesting parallels between systems of interacting non-Abelian anyons and quantum mag- 
netism, occuring in ordinary SU(2) quantum magnets. Here we consider theories of so-called su(2)fc anyons, 
well-known deformations of SU(2), in which only the first k + 1 angular momenta of SU(2) occur. In this 
manuscript, we discuss in particular anyonic generalizations of ordinary SU(2) spin chains with an emphasis 
on anyonic spin S = 1 chains. We find that the overall phase diagrams for these anyonic spin-1 chains closely 
mirror the phase diagram of the ordinary bilinear-biquadratic spin- 1 chain including anyonic generalizations of 
the Haldane phase, the AKLT construction, and supersymmetric quantum critical points. A novel feature of 
the anyonic spin-1 chains is an additional topological symmetry that protects the gapless phases. Distinctions 
further arise in the form of an even/odd effect in the deformation parameter k when considering su(2)k anyonic 
theories with k > 5, as well as for the special case of the su(2)4 theory for which the spin-1 representation plays 
a special role. We also address anyonic generalizations of spin- 1/2 chains with a focus on the topological pro- 
tection provided for their gapless ground states. Finally, we put our results into context of earlier generalizations 
of SU(2) quantum spin chains, in particular so-called (fused) Temperley-Lieb chains. 

PACS numbers: 05.30.Pr, 03.65.Vf, 03.67.Lx 
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I. INTRODUCTION 

Ever since the early days of condensed matter physics, 
quantum magnets have played an integral role in shaping our 
understanding of interacting quantum many-body systems. 
Following the experimental discovery of the high-temperature 
superconductors whose undoped parent compounds typically 
are antiferromagnets, the study of quantum magnets has fur- 
ther intensified yielding a plethora of deeper insights. Early 
on, quantum spin chains - typically one-dimensional arrange- 
ments of SU(2) spins - have become prototypical systems 
that proved to be fruitful ground for analytical descriptions 
and quasi-exact numerical analysis^. One seminal result was 
the exact solution of the antiferromagnetic spin- 1/2 Heisen- 
berg chain via the Bethe ansatz and its description in terms of 
conformal field theory. Another crucial contribution was Hal- 
dane's realization 2 that the antiferromagnetic spin-1 Heisen- 
berg chain forms a gapped state with characteristic zero- 
energy edge states for open boundary conditions - a princi- 
ple observation that holds true for all half-integer and inte- 
ger spin chains. More recently, it has been found that the 
gapped Haldane phase of the spin-1 chain is an example of 
a symmetry protected topological phase 34 making it a one- 
dimensional cousin of topological insulator states in two and 
three dimensions 5 , which have attracted much recent interest. 

Over the years, a plethora of physical systems that con- 
nect to the elementary physics of quantum spin chains have 
been identified, including transition metal oxides 6 , Au quan- 
tum wires on semiconducting surfaces^, or ultra-cold atoms 
in optical lattices 8 . Recently, it has been realized that cer- 
tain 'deformations' of quantum spins can be used to describe 
some of the more peculiar topological properties of exotic 
quasiparticles, so-called non-Abelian anyons, that arise in 
certain topologically ordered systems, including certain frac- 
tional quantum Hall states 9 , p x + ip y superconductors^, het- 
erostructures of topological insulators and superconductors 11 , 
heterostructures of spin-orbit coupled semiconductors and 
superconductors 12 and possibly certain Iridates 13 which may 
effectively realize the Kitaev honeycomb model 14 . To be more 
specific, the deformations of quantum spins are representa- 
tions of the anyon theories called su(2)fc, which can be de- 



scribed as theory of ordinary SU(2) quantum spins that is de- 
formed in such a way that only the first k + 1 (generalized) 
angular momenta 

1 k 

j = o, 2 ,i,...,- 

can occur. These generalized angular momenta capture the 
non-Abelian properties of the anyonic quasiparticles present 
in the su(2)^ theory. For instance, the non-Abelian nature of 
the so-called Majorana fermion is captured by the general- 
ized angular momentum 1/2 of the su(2)2 theory. The same 
holds for so-called Ising anyons, while Fibonacci anyons can 
be represented by the generalized angular momentum 1 of the 
su(2) 3 theory. Similar to the coupling of two ordinary spins, 
a pair of generalized angular momenta can be combined (or 
'fused') into a new set of joint quantum numbers. For in- 
stance, for k > 2, two generalized angular momenta 1/2 can 
be combined to form either a state with generalized angular 
momentum or a state with generalized angular momentum 
1, which is written as 

1/2x1/2 = + 1, (1) 

reminiscent of two ordinary spin 1/2's forming either a singlet 
or triplet state. Similarly, two generalized angular momenta 1 
can be combined into 

1x1 = + 1 + 2 (2) 

for deformation parameters k > 4. For lower values of k, the 
rules differ because the number of representations is limited 
by k. In particular, for k = 3, one finds 1x1 = + 1, 
while for k = 2, one has 1x1 = 0. Finally, for k = 1, the 
general momentum 1 is not allowed. For the anyonic theories, 
the above equations are often referred to as fusion rules. 

The many-body physics of a set of interacting non-Abelian 
anyons can be captured by a Hamiltonian that is formed by 
pairwise interactions which assign energies to the different 
outcomes in the above fusion rules. Such an approach is a 
straightforward generalization of the conventional Heisenberg 
model, whose pairwise interaction term JSi ■ Sj is simply a 
projector onto the singlet state, which is energetically favored 
for antiferromagnetic couplings (J < 0) or penalized for fer- 
romagnetic couplings (J > 0). 

The first step in this direction was taken by some of us for 
anyonic spin- 1/2 chains in Ref. (PTBT l and later generalized to 
spin-1 chains in Ref. ( 17) by the current group of authors. The 
careful analysis of the ground states of these one-dimensional 
systems has resulted in a number of insights. First, anyonic 
spin-1/2 chains typically form gapless ground states which 
can be described in terms of conformal field theory 16 . These 
gapless states turn out to be protected by a topological sym- 
metry inherent to the anyon chains that renders them stable 
against local perturbations 16 18 . Moreover, these gapless states 
can in fact be interpreted as edge states that reveal the true 
ground state of a two-dimensional set of anyons - a novel 
topological liquid that is separated by the original topologi- 
cal liquids (of which the anyons are excitations) by an edg&^l 



3 



This picture has been verified by a careful analysis of ladder 
systems, in which multiple chains are coupled-^. 

Going beyond spin- 1/2 chains, we began to study the 
physics of anyonic spin-1 chains with first results be ing re- 
ported in a preceding (much more condensed) papep^l In 
the manuscript at hand, we provide an in-depth discussion of 
these anyonic spin-1 chains. We find that many of the distinc- 
tive features of ordinary SU(2) spin- 1/2 and spin-1 chains also 
hold for their anyonic cousins. For instance, anyonic spin-1 
chains exhibit a gapped topological phase for antiferromag- 
netic couplings - the anyonic generalization of the Haldane 
phase. Exploring the phase diagram of chains of pairwise 
interacting spin-1 anyons, we find a striking resemblance of 
the anyonic phase diagram to the one of the ordinary bilinear- 
biquadratic spin-1 chain. In particular, we find multiple gap- 
less phases (and phase transitions) in addition to the gapped 
Haldane phase. For the former, a similar topological protec- 
tion mechanism and edge state interpretation holds as for the 
gapless phases of the anyonic spin- 1/2 chains^. 

The focus of this manuscript is to provide an exhaustive de- 
scription of the phase diagram(s) of the anyonic spin-1 chains. 
Our exploration of these systems has led to a large amount of 
results as the phase diagrams turned out to be much richer 
than initially anticipated. In particular, we find two families 
of phase diagrams depending on whether the deformation pa- 
rameter k of the su(2)fc anyonic theories is odd or even. More- 
over, we obtain a distinct phase diagram for k — 4, a result 
that can be explained by the special role played by the gener- 
alized angular momentum 1 in the su(2)4 theory. 

In order to guide the reader through these various results 
we have taken some care to structure the manuscript as fol- 
lows: We will start with an introduction to the anyonic su(2)fe 
theories and a description of the anyonic generalization of the 
Heisenberg model in Section II. The following sections will 
then give a detailed expose of our results, devoting Sec. Ill 
to the discussion of anyonic spin-1 chains with odd deforma- 
tion parameters k > 5, followed by a discussion of the case of 
even deformation parameters k > 6 in Sec. IV. In Sec. V we 
will turn to the case of k = 4 for which the spin-1 representa- 
tion plays a special role and a rich phase diagram is obtained. 
We will then turn to anyonic spin- 1/2 chains and discuss their 
physics, in particular their topological stability in Sec. VI. We 
will end with a broader discussion of our results, in particu- 
lar in light of other deformations of conventional spin chains 
such as continuous su(2) 9 deformations or so-called (fused) 
Temperley-Lieb spin chains. The main part of the manuscript 
is followed by an appendix that provides the technical details 
of our calculations. 



II. THE ANYONIC QUANTUM SPIN CHAIN 
HAMILTONIANS 

In light of the recent interest in topological phases of matter, 
it is of great importance to gain an understanding topological 
models in their simplest incarnation, and we will thus study 
one-dimensional chains of interacting non-Abelian anyons. In 
this section, we will briefly explain the models by drawing 



parallels with ordinary one-dimensional spin chains. More- 
over, we will explain why the 'topological' nature of these 
models goes beyond the fact that they are constructed from 
'topological' particles, namely non-Abelian anyons. 

One of the prototypical one-dimensional spin chain models 
is the Heisenberg model, in which SU(2) spins interact via a 
'spin-spin' interaction of the type S{ ■ Sj, where the labels i 
and j denote the locations of the interacting spins. Often, one 
restricts the interaction to nearest-neighbor, or next-nearest- 
neighbor pairs of spins. For the description of the anyonic 
quantum spin chains, it will be beneficial to think of this in- 
teraction in terms of the total spin of the two interacting spins. 
In this paper, we will only consider nearest neighbor interac- 
tions. 

As a first example, we look at conventional SU(2) spin- 1/2, 
and consider the total spin St — (Si+Si+i) of two interacting 
spins Si and Si+i, whose magnitude is characterized by the 
eigenvalue of (St) 2 = (Si + Si+i) 2 . Because the total spin 
St can be either or 1, with S T eigenvalues and 2, we can 
write 

(5 l + 5 l+1 ) 2 = 0P l (0) +2^ (1) , (3) 

where the projection operator P> ' projects onto the total spin 
s channel of the two spins Si and Sj+i. Evaluating the left 
hand side, one obtains 

& ■ S i+1 = if> - 3 -h = -P^ + \h , (4) 

where in the last step we used that we can rewrite the iden- 
tity operator as Ij = P^ + P^\ which holds in the case of 
spin-1 /2. We conclude that the Heisenberg interaction assigns 
energy to two interacting spins, depending on their combined 
spin, and the Heisenberg Hamiltonian can be written in terms 
of projectors as 

H = JJ2P{°\ (5) 

i 

where J = 1 corresponds to an antiferromagnetic coupling, 
and J = — 1 to the ferromagnetic version. 

For spin-1, one can similarly write the bilinear and bi- 
quadratic terms Si-Si + i and (Si-Si + i) 2 respectively, in terms 
of the projection operators and P^ . In particular, the re- 
lations 

(S t + S t+1 f = 2P l {1) +6Pf } 
(S ? , + S ? m ) 4 -4^ (1) +36^ (2) , (6) 
can be rewritten as 

(5,-5 l+1 ) = P l (1) +3^ (2) -2I l 
(Si ■ S t+1 ) 2 = -3P 2 (1) - 3Pf } + 41, . (7) 
Consequently, the bilinear-biquadratic spin-1 Hamiltonian 

H hh = c o s (0bb)(^ • S i+1 ) + sm(e hh )(§i ■ S l+1 ) 2 (8) 
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can be expressed in terms of the projectors pj 1 ' and P^ as 
follows, 



i 

= ^cos 9 2A P l (2) -sin0 2 ,i P i 



(i) 



(9) 



Here, the relation between the two angles 9 2 i and 6*bb is given 
by 



tan6> bb -l/3 tan0 2 ,i + l/3 

tan U21 — t. tan ^bb = ' t, 

1 — tan fbb 1 + tan 6*2.1 



(10) 



We will now shift our attention to anyonic degrees of free- 
dom. Details about anyon models, in particular those of type 
su(2)fc, can be found in appendix [A] A general introduction 
can be found, e.g., in references Ref. (114115120b . Here, we 
will only introduce those concepts that are necessary for defin- 
ing the chain Hamiltonians. The Hamiltonians for the anyon 
chains that we will consider in this paper are of the form of 
Eq. |9j. The projectors P^ in that equation have however a 
different meaning for anyons (as compared to ordinary spins) 
which will be defined in Eq. ( fT2| i below. 

Anyons are labeled by generalized angular momenta, or - in 
the language of anyons models - 'topological charges'. These 
generalized angular momenta correspond to quantum num- 
bers, just as in the case of ordinary spin degrees of freedom. 
The notion of combined spin, or tensor product of spins, cor- 
responds to the notion of 'fusion' in the language of anyons, 
and can in general result in more than one type of anyon. The 
possible outcomes are called 'fusion channels'. The general- 
ization of the Heisenberg interaction for spins to the anyonic 
case is to assign an energy to two interacting anyons based 
on their fusion channel. How this is done in practice, will be 
described in more detail below and in appendix [C] 

The class of anyons considered in this paper is derived from 
SU(2) where spin-S ranges from 5" = 0, 1/2, 1,3/2, .... In 
contrast, su(2)fc anyons contain only a subset of generalized 
angular momenta, namely 

1 k 

j = o, 2 ,i,...,-. 

The truncation, characterized by the 'level' k, has two impor- 
tant consequences which we will describe in the following. 

The first consequence concerns the fusion rules of the 
anyons. The tensor product of two SU(2) spins Si and S 2 
decomposes as 



associative. In particular, the fusion rules of su(2)fc anyons 
read 

J1XJ2 = I Ji — J2 1 + (lii — is i +1) H \-mm(j 1 +j 2 ,k-ji-j 2 ). 

(11) 

The second important consequence of the truncation fol- 
lows from the fusion rules. The dimension of the Hilbert 
space of a number A" of ordinary SU(2) spin-l/2's is equal 
to 2 N , and the spins can add up to a maximum spin of N/2. 
In contrast, the dimension of the Hilbert space of a number A" 
of j = 1/2 anyons in the su(2)^. theory is smaller than 2 N . 
In appendix [A] it is shown that the dimension of the Hilbert 
space for N j — 1/2 anyons grows as dy 2 , asymptotically for 

large N, where di/2 — 2 cos fpp^ is the so-called quantum 

dimension of the j = 1/2 anyon. For 1 < k < 00, this im- 
plies that the effective number of degrees of freedom for each 
anyon is irrational. This is less mysterious than it sounds: all 
this is saying is that one can not think of the Hilbert space of 
A" anyons as a tensor product of A" one-anyon Hilbert spaces. 

Because the Hilbert space does not have a tensor product 
structure, an alternative description of the state space and the 
Hamiltonian acting on it is needed. We will describe here 
how this can be done, but leave the details for the appendices 
where we also give an explicit description of the Hamiltonians 
studied in this paper. 

The Hilbert space of a chain of anyons can be described in 
terms of a so-called 'fusion tree'. In Figure [T] the fusion tree 
for a chain of 'spin-1' anyons is displayed. The lines in the 
fusion tree carry a label indicating the type of anyon the line 
corresponds to. The lines coming from above correspond to 
the spin-1 anyons which constitute the chain. The horizontal 
lines, labeled by Xi, are the actual degrees of freedom. The 
possible 'values' of the Xi are the same as those of the anyons 
present in the anyon model, namely Xi = 0, 1/2, . . . , k/2, in 
case of su(2)fc anyons. The Xi cannot be chosen arbitrarily, 
but may only take values such that the fusion rules are obeyed 
at the trivalent points. For example, the anyon type xi has to 
appear in the fusion product of xq x 1, the anyon type x 2 ap- 
pears in the fusion product x\ x 1, and so on. Each labeling 
of the fusion tree that is consistent with the fusion rules corre- 
sponds to a (orthonormal) state in the Hilbert space, and these 
states span this space. 

Typically, we will use periodic boundary conditions xl = 
xq, which implies that xq has to appear in the fusion prod- 
uct x 1, where L denotes the number of sites of 
the chain. States in the Hilbert space will be written as 
\x ,xi, . . .,x L -i). 

The Hamiltonian assigns an energy based on the fusion 
channel of two neighboring anyons in the chain. However, 



Si ® S 2 = I Si - s 2 



© (Si + S 2 ) 
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The process of taking tensor products is associative, and the 
same is true for the fusion rules. Because of the truncation 
in the su(2)fc theory, the SU(2) tensor product rule has to be 
modified. It turns out that there is only one way of doing 
this, consistent with the requirement that the fusion rules are 



X Xi x 2 x 3 ... 
FIG. 1: The anyonic spin-1 chain. 
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| N f^2 | | 
Xq Xl X2 X3 . . . Xq X\ X3 . . . 

FIG. 2: The basis transformation for the anyonic spin-1 chain. 



in the above discussed representation of the Hilbert space (see 
Figure [TJ, the fusion channel of two neighboring anyons is 
not explicit. To remedy this problem, we employ a local basis 
transformation which changes the order in which the anyons 
are fused. This is permissible because of the associativity of 
the fusion rules. For ordinary SU(2) spins, this basis transfor- 
mation is described in terms of the Wigner 6j-symbols. In the 
case of anyons, this basis transformation is described by what 
are known as the F-symbols. A detailed discussion of the F 
symbols, as well as explicit representations for su(2)j. anyons 
can be found in appendix [B] 

The basis transformation is depicted in Figure [2] On the 
left hand side, x\ is fused with a spin-1 anyon, resulting in 
X2, which is subsequently fused with the next spin-1 anyon, 
resulting in anyon type x s . After the basis transformation, 
one first fuses the two spin-1 anyons, resulting in x%, which is 
fused with x\, resulting in the anyon type 2:3. Both bases are 
equivalent; however, in the second basis, the fusion channel 
of the two spin-1 anyons is explicit, namely x 2 . Thus, after 
performing this basis transformation, one can assign the ap- 
propriate energy based on the value of X2- Subsequently, one 
transforms back to the original basis. The operator projecting 
onto the anyon— j channel of two neighboring anyons i and 
i + 1 is thus given by 



0) : 



(12) 



where F, is shorthand for the local basis transformation de- 
picted in Figure^ The operator 11^ projects onto the fusion 
channel x^ = j, i.e., the fusion of two anyons into an anyon 
of type j is penalized with energy E = 1, while the other pos- 
sible fusion channels are assigned E = 0. For explicit matrix 
representations of P^ we refer to appendix |c| 

It is important to realize that the form of the projector ( fT~2| > 
is universal and applicable to anyonic chains composed of ar- 
bitrary types of anyons. Changing to a different anyon model 
will merely result in a different structure of the Hilbert space 
and different F-symbols. 



A. Topological symmetry 

In this section, we present a detailed discussion of the 
'topological symmetry operator'. The Hamiltonians consid- 
ered in this paper commute with the topological symmetry 
operator, and the associated symmetry plays a crucial role in 
the analysis of the anyonic chain models. In 'equation' ( [13) , 
a chain of type-j anyons with periodic boundary conditions is 



displayed (in this particular case, L = 3). 



•''0 




j 




J 






Xl 



Xq 




(13) 



For each type of anyon I, there exists a topological operator Y\. 
The action of this operator Yj on the state |xrj, £1, • • • , xl-i), 
displayed in ( fT3j ) for L = 3, can be described as follows. First, 
an additional anyon of type / is created inside the spine of the 
fusion tree, as displayed in ( |14) , 



Xn 
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J 




X2 


c ^ 



Xq 




(14) 



This additional spin-^ anyon is 'merged' with the fusion dia- 



gram by applying an F-matrix, namely, ( F l 



7 x ,x ,l\ 



I 



ing in the state 

j2( F r xo ' 1 )^ 1*0,31,..., *i-i> 

x Q 

as depicted in (TT3). 



, result- 
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Xl 
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(15) 



Next, one 'moves' the additional spin-Z anyon around the ring, 
by applying additional F-matrices. After the first step, one 
obtains the state 

{Fir l )S , xl ,..., XL .i) 



as illustrated in (116 




(16) 



Another move of this sort gives 



XL-l) 
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as shown in ( fTT| i. 




(17) 



Finally, after L steps, one has come full circle, giving rise to 
the states 



L-l 



, e (^"" j )r n (*r + "');r 

as depicted in ([18^, for L = 3. 



X ,Xi, . . . ,X L -l) 





3 


3 


3 






x 2 


r" 











r 





(18) 



From the general properties of anyon models (see e.g. Ref. ( 
114115b ). we find that x'q = x' Q (the overall topological quan- 
tum number of an isolated set of anyons can not change). We 
can now remove the additional spin-/ anyon in the same way 
as we added it, thereby finishing the operation of acting with 
Yi on the state \xo, x\, . . . ,xl-i)- Thus, we obtain the ex- 
pression 

Yi \x ,xi, . . . ,x L -x) = 

E II (/••;•' ' ■')'/' -'w, <•/. , . d9) 



We can now state the matrix elements of the topological oper- 
ator Yi in the fusion tree basis 



L-l 



(x' ,x[, . . . ,x' L _j_\ Yi \x ,Xi, . . .,X L -!} 



n 



i=0 



Xi + l 

(20) 

The above definition of the topological operator does not 
depend on whether the additional spin-/ anyon is encircled by 
the anyon chain (as in Figures (fT4|i-(fl"8]l) or whether the ad- 
ditional spin-/ anyon encircles the entire anyon chain. When 
using the latter description of the topological operator, one can 
think of the additional spin-/ anyon as going around the 'fu- 
sion product' of all the spin-j anyons constituting the anyonic 
chain, or better, encircling the flux through the chain. This 
flux through the chain is related to the additional spin-/ anyon 
as follows, 



(21) 




where i denotes the flux going through the chain, and the ma- 
trix S is the modular S'-matrix of the anyon model. For a 
derivation of Eq. pTj ), see e.g. Ref. J14I15I) . and the explicit 
form of S in the case of su(2)fc anyons is given in appendixfA] 
The definition of the topological operator contains elements 
of the F matrices only. This is also true for the anyonic spin 
Hamiltonians we consider in this paper. It follows that the 
operators Y\ commute with the Hamiltonian and that a topo- 
logical quantum number can be assigned to all the eigenstates. 
This has far reaching consequences for the stability of the crit- 
ical phases. Excited states which are relevant in the renormal- 
ization group sense (i.e., have energy smaller than 2) may lie 
in a different topological sector than the ground state and thus 
do not drive the system into a different phase. In addition, 
we will see that the operators Yi play an important role in the 
zero-energy ground states at the AKLT point in the Haldane- 
gapped phase of the spin-1 models. 



III. ANYONIC SU(2) fc SPIN-1 CHAINS: ODD k > 5 
A. Introduction 

We will start our discussion of anyonic quantum spin chains 
with the anyonic version of the ordinary SU(2) spin-1 Heisen- 
berg chain, which has long been appreciated as one of the 
paradigmatic spin chain models. For antiferromagnetic cou- 
plings the spin-1 chain is well known to form a gapped phase, 
in distinction from the gapless spin- 1/2 Heisenberg chairP. 

In the following sections, we discuss in detail the anyonic 
su(2)/ £ deformations of the ordinary SU(2) spin-1 chain. We 
will see that much of the seminal features of the SU(2) spin-1 
chain carry over to these anyonic deformations with a number 
of new subtleties arising. One is a dependence of the observed 
phases and phase diagrams on the deformation parameter k. 
In particular, we find an even/odd effect in k for k > 5 and 
a distinctive behavior for k = 4. We have therefore split our 
discussion of the anyonic spin-1 chains into three different 
sections. We will address anyonic spin-1 chains with odd k > 
5 in the remainder of this section, in which we will also give 
a brief recount of the phase diagram of the ordinary SU(2) 
spin-1 chain. The subsequent section will be devoted to the 
case of k > 6 with k being even. Finally, an entire section 
is devoted to a detailed account of the physics for the special 
case of k = 4. 



B. The ordinary SU(2) Heisenberg spin-1 chain 

Before addressing the physics of the anyonic spin chains we 
briefly recapitulate the phase diagram of the ordinary SU(2) 
spin-1 Heisenberg chain. While the latter is typically dis- 
cussed as a circle phase diagram in terms of bilinear and bi- 
quadratic spin exchange, we will recast the phase diagram in 
terms of the projector representation in Eq. |9]) - the generic 
representation of anyonic spin chains. Fig.[3]shows the phase 
diagram in the projector representation of Eq. |9]). It contains 
four different phases, of which two are gapped phases and two 
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SU(2) spin-1 chains 

SU(3) 




SU(3) 



FIG. 3: (color online) Phase diagrams of the ordinary SU(2) spin-1 
chain in a projector representation (|9j with J\ — — sin(#2,i) and 
J2 = cos(02,i). 



are gapless phases. The well known Haldane phase 2 extends 
in the parameter regime — arctan(2/3) < 6*2.1 < tt/2 and 
includes the so-called Affleck-Kennedy-Lieb-Tasaki (AKLT) 
poinP^at #2,1 = (in which only the projector is present 
in the Hamiltonian), at which the exact form of the ground- 
state wave function in terms of a valence bond solid state 
can be obtained. The conventional (gapped) Heisenberg chain 
(bilinear in spin-1 operators) with antiferromagnetic coupling 
corresponds to 2 .i = — arctan(l/3). The second gapped 
phase is a (spontaneously) dimerized phased that occurs in 
the parameter regime — tt/2 < 6*2,1 < — arctan(2/3). The 
phase transition at #2,1 = — arctan(2/3) between the two 
gapped phases is described by the su(2) 2 conformal field the- 
ory with central charge c = 3/2, which happens to possess 
N = 1 supersymmetry - a result that can be obtained by 
means of a (nested) Bethe Ansatz^l. 

At the other end of the Haldane gapped phase, 6*2,1 = 
tt/2, there is a phase transition to gapless phase that ex- 
tends over the range tt/2 < 6*2.1 < 37r/4. This critical 
phase can be described by a conformal field theory with cen- 
tral charge c = 2. There are characteristic quadrupolar (ne- 
matic) spin correlations 24 in this phase, as well as a three sub- 
lattice structure 2 ^ resulting in soft modes at momenta K = 
0, 27r/3, 47r/3. At the transition from the gapped Haldane 
phase to this critical nematic phase at #2,1 = tt/2, the system 
has enhanced SU(3) symmetry. This point in the phase dia- 
gram of the spin-1 SU(2) chain represents actually the SU(3) 
chain with a fundamental representation at each site, which is 
known to be described by the SU(3)i conformal field theory. 
(This chain is again exactly solvable by a Bethe Ansatz 2627 .) 

Finally, there is a gapless ferromagnetic phase, extending 
over the parameter range 37r/4 < #2,1 < 37r/2. The phase 
transitions from this phase to both the adjacent dimerized 
phase as well as the nematic phase are first order. In the vicin- 



su(2)k spin-1 chains (odd k) 




FIG. 4: (color online) Phase diagrams of the anyonic su(2)^ spin- 
1 chain with odd A; in a projector representation ^ where Ji = 
— sin(#2,i) and J2 = cos(#2,i). With increasing (odd) index k > 5 
the phase boundaries move as indicated by the arrows. 



ity of the transition between the dimerized and ferromagnetic 
phase, early analytical work 28 suggested the possibility of an 
intermediate nematic p hase, w hich, however has later been 
found to not materializ e 1 24 ! 29 ! 30 !. 



C. Phase diagram of the anyonic spin-1 chains - overview 

In this section, we provide an overview of the phase dia- 
gram of the anyonic spin-1 chains for odd k > 5. This phase 
diagram bears great resemblance to the corresponding phase 
diagram of the SU(2) spin-1 Heisenberg chain (Fig. [3}. The 
generic phase diagram for the su(2)^ spin-1 chain is given 
in Figure [4] In Figures [5] and [6] we display the phase di- 
agrams for k — 5 and k = 7, as well as the character- 
istic spectra of the four different phases and the (N = 1) 
super-symmetric critical point which separates the Haldane 
gapped phase and the phase which will be called "Z 2 sublat- 
tice phase" (this is the phase intervening between the Haldane 
phase and the Z^-parafermion phase, and it encompasses the 
angles 2 ,i < -0.19tt w - arctan(2/3)). 

The spin-1 anyonic spin chain is gapped in a finite region 
around 6*2,1 = 0. This gapped phase is the anyonic analogue 
of the Haldane gapped phase, and the point #2,1 = is equiv- 
alent to the AKLT point. At this point, the Hamiltonian penal- 
izes the fusion of two neighboring anyons in the spin-2 chan- 
nel. The ground states with periodic boundary conditions can 
be found exactly at this point, for all k, and the ground state 
degeneracy is (k + l)/2. 

For 02,i < 0, there is a phase transition at 6*2,1 ~ — 0.197T 
into an extended critical region. The position of this phase 
transition did not show any appreciable dependence on the 
value of k (remember that k > 5 throughout this section). 
This gapless region occurs where the ordinary SU(2) spin-1 
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chain is in the gapped dimerized phase. This difference in be- 
havior is the most remarkable distinction between the ordinary 
SU(2) spin-1 chain, and the anyonic spin-1 chains. 

The critical point at 2 ,i ~ — 0.197T w — arctan(2/3), sep- 
arating the Haldane phase and the extended critical region, is 
described in terms of an N = 1 super-symmetric minimal 
conformal model. 

For angles 02. 1 > 0, there is a phase transition from the Hal- 
dane phase into another extended critical region which bears 
some resemblance to the extended nematic region in case of 
the ordinary spin-1 chain. In particular, this phase has a Z3 
sublattice structure. The location of the phase transition does 
depend on k, and moves towards 82 1 = tt/2 with increasing 
k. 

Finally, there is an extended critical region in the vicinitiy 
of 02,i = ?r> the point where the fusion of two neighboring 
anyons into the spin-2 channel is favored. This critical phase 
is the anyonic analogue of the ferromagnetic phase of the or- 
dinary spin-1 chain, and the critical behavior is described by 
the Zk parafermion conformal field theory. 

The phase transitions from the ferromagnetic phase to the 
neighboring extended critical regions are first order. The 
phase transition into the anyonic version of the nematic phase 
occurs at 02.1 = 37r/4, independent of the value of k. The lo- 
cation of the other phase transition depends on k, and moves 
towards 02, 1 = 37r/2 for increasing k. 

Below, we will discuss in detail each of the phases men- 
tioned above. We will focus on the topological properties and 
the similarities to the ordinary SU(2) spin-1 chain. 



D. Critical phases 

We investigate the phase diagram of our model numeri- 
cally using exact diagonalization. In our analysis, we follow 
a standard procedure to determine the conformal field theory 
describing the behavior of the extended critical regions and 
the critical points: the numerically obtained spectrum is first 
shifted (by some constant offset) such that the ground state 
has zero energy zero. The spectrum is then rescaled such that 
the energy of the lowest lying excitation matches the energy 
of the lowest lying excitation of the conformal field theory de- 
scribing the phase. The so obtained energy spectrum is finally 
compared to the energy spectra of candidate CFTs. The CFT 
(if any) which matches the numerically obtained energy levels 
is the one describing the system at the angle 0. We note that 
the list of candidate CFTs is limited: If the chain is critical, 
each energy level in the spectrum corresponds to a field in the 
applicable conformal field theory. These fields satisfy fusion 
rules which have to be compatible with the fusion rules of the 
underlying su(2)^, theory. This constraint restricts the candi- 
date conformal field theories that could describe the criticality 
of anyonic quantum chains. 

The eigenenergies in a system of finite size described by a 
conformal field theory take the forrrP2 



where the velocity v is an overall scale factor, and c is the cen- 
tral charge of the CFT. The scaling dimensions h + h take the 
form h = h° +n,h — h° + n, with n and n non-negative inte- 
gers, and h° and h° are the holomorphic and antiholomorphic 
conformal weights of the primary fields in the given CFT. The 
momenta K (in units 2ir / L) are such that K = h—h + Kq or 
K = h — h + Kq + L/2, where K is a constant shift of the 
momentum that determines at which momentum the primary 
field occurs. This shift can be determined from the numer- 
ics, and is not fixed by conformal symmetry. Thus, different 
microscopic realizations of the same conformal field theories 
can give rise to different values for Kq. 



As explained in section II A the anyonic spin chains have a 
topological symmetry; all the states in the spectrum can there- 
fore be assigned a topological quantum number. The pos- 
sible eigenvalues of the topological symmetry operator, also 
denoted as topological quantum numbers, are in one-to-one 
correspondence with the types of anyons which appear in the 
particular anyon theory considered. 



1. Zk-parafermion phase 

We begin the discussion of the phase diagram given in Fig- 
ure[4]with the Zk -parafermion phase which corresponds to the 
gapless ferromagnetic phase in the SU(2) spin-1 chain. In the 
anyonic spin-1 chains, this phase contains the point 2 .i = n 
where it is favorable for two neighboring anyons to be in the 
spin-2 channel. One of the phase boundaries of this phase 
is located at 2 1 = 37r/4. The location of the other phase 
boundary depends on k: with increasing k, it moves towards 
the location of the phase boundary in the SU(2) spin-1 chain 
(at angle 9 2 ,i = 3tt/2). 

The spectra at angle 2 i = ir for k = 5 and k = 7 are 
displayed in the middle panel of Figures [5] and [6] respec- 
tively. The energy spectra were rescaled such that the energy 
of the lowest excitation matches the energy predicted by the 
Zk parafermion conformal field theory^. Some details of this 



CFT are reviewed in appendix E4 In the Figures, we indi- 
cate the locations of the energies of the states corresponding 
to the primary fields by green squares, while blue crosses cor- 
respond the numerically obtained energy levels. We find good 
agreement between numerically obtained energy spectra and 
the Zk parafermion CFTs for both the su(2)s and su(2)7 anyon 
models. For su(2)5, we also indicate the location of a few 
descendant fields that match the numerical prediction. Gen- 
erally, the identification of descendant fields is more difficult 
due to finite size effects. 

The fields of the Zk parafermion theory carry two la- 
bels, (l,m) that take the values / = 0, 1, ... (fc — 1), and 
m = 0, 2, . . . 2(k — 1). The momentum and topological quan- 
tum number of the fields is determined by the labels m and 
I, respectively. The topological quantum number simply is 
given by I. For the momentum, the following relation holds: 

is 2mir 

fe ' 
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FIG. 5: (color online) The su(2)s spin-1 chain: The energy spectra for the various phases of the phase diagram are shown in the upper left 
panel. For the critical phases/point the energy spectra have been rescaled to match the conformal field theory prediction given in Eq. \22\ . 
Green squares indicate the location of the primary fields, red circles the descendant fields. The energies predicted by conformal field theory 
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We find that there are no relevant primary fields which have 
the same set of quantum numbers as the identity field. This 
implies that there are no relevant operators that can be added 
to the Hamiltonian to drive a phase transition if both transla- 
tional and topological symmetry are left unbroken. This phase 
is an example of a critical phase whose criticality is protected 
by the topological symmetry. 



indicate the topological sectors of some of the low lying fields 
and give the scaling dimensions of the primary fields. 



4. Superconformal critical point 



2. Zi-phase: (A, D) modular invariant of coset 
su(2) k -x X su(2)i/su(2) k 

Upon increasing 6*2,1, one encounters a first order transition 
from the Zk parafermion phase into a different extended criti- 
cal phase that has a Z 2 sublattice symmetry. We identified the 
CFTs describing these critical phases for k = 5 and k = 7 
as Virasoro conformal minimal models, with central charge 
c = 1 — ( fc+1 ) 6 ( fc+2 ) • However, the field content describing 
the criticality is not the 'usual' minimal model - the diagonal 
(A, A) modular invariant - but the so-called (A, D) modular 
invariant which contains a different number of fields. Details 
of these different modular invariants can be found irPUEll F r 
our purposes, it suffices to notice that some of the primary 
fields in the (A, A) invariant do not appear in the (A, D) in- 
variant while others appear twice. The details of this CFT are 
summarized in table VII in appendix |E 1| T he scaling dimen- 
sions of the fields are given in equation ( |E1[ >. 

Again, it is possible to identify the topological sectors and 
the momenta at which the various fields occur from the labels 



of the fields. As discussed in appendix E 1 the fields can be 
labeled by (r, s), where s takes the values s = 1,3, ... ,k. 
The topological sector is determined by (s — l)/2, while the 
momenta are fixed by the r label. In particular, for k = 5, the 
fields with labels r = 1, 5 occur at K = 0, while the fields 
with label r = 3 occur both at K = 0,ir. For k = 7, the fields 
with r = 1, 3, 5, 7 occur at K = 0, while the fields at r = 4 
are doubly degenerate and occur K = tt. 



3. Z^-phase: coset su(2)k-4 x su{2)4,/ su(2)k 



The transition between the Z2 phase and the Haldane 
gapped phase occurs at the angle 6*2.1 ~ — 0.197T, which 
shows little dependence on the level k. The critical point itself 
is described by a N = 1 superconformal minimal model 38 , 
su(2)2 xsjj(2 ) fc _2 ^ rj e t a ii s on m is theory can be found in ap- 



pendix E 2 In the limit of k — > 00, this theory approaches the 



su(2)2 theory, which describes the critical point in the SU(2) 
spin-1 bilinear-biquadratic spin chain. 

In the spectra for k = 5 and k = 7 of the anyonic spin- 
1 chain at this critical point, we indicate the scaling dimen- 
sions and topological sectors of the primary fields which are 
labeled by (r, s). Like in the other coset models (excluding 
the Zk parafermion theory), the label s is associated with the 
su(2)fc denominator of the coset and hence labels the topolog- 
ical sector. The momentum at which the primary fields appear 
is determined by K = (r + s mod 2)7r. 

The superconformal critical point separates the Haldane 
gapped phase from the Z2 sublattice critical region. There- 
fore, we expect that there will be a relevant perturbation which 
drives the phase transition between these two different phases, 
and that this perturbation does not break any symmetries. A 
relevant perturbation is a field which has the same quantum 
numbers as the ground state and whose scaling dimension is 
smaller than two. Such a field indeed exists: it carries the la- 
bels (r, s) = (3, 1) and has scaling dimension 1 + $ , i.e., it 
is a relevant field for all fc. We note that at K = tt, there also 
is a relevant field with labels (r, s) = (2, 1) which has scaling 
dimension | + Jj. As a consequence, a gap is expected to 
develop if a perturbation which staggers the chain is added to 
the system. 



At 6*2.1 = 37r/4, there is a first order transition between 
the Zk 'ferromagnetic' phase and a critical region of critical- 
ity that exhibits a Z 3 sublattice symmetry. We determined 
that the CFT describing the Z3 critical region is a series of 

sii(2) 4 xs«(2) fc _ 
su(2) k 



coset models with S3 symmetry, namely 3 " w t^/o" Wfc 4 . In 



appendix E3 we list some details of these coset models, in 
particular, the scaling dimensions of the primary fields (a de- 
tailed analysis can be found i n 36 * 37 ! ). The primary fields are 
labeled by two integers (r, s). As was the case for the Z2 crit- 
ical phase, only a subset of the fields appear in the spectrum, 
namely those with r + s even. In addition, the label s has to 
be odd, and it determines the topological quantum number via 
(s — 1) /2. The location of the second endpoint of this Z3 crit- 
ical region (i.e., the transition to the Haldane gapped phase) is 
found to vary with k. 

In Figures [5] and [6] we display representative energy spec- 
tra for this phase (angle 6*2,1 = 0-7%). In these spectra, we 



J. Stability of the critical phases 

We recapitulate that in all three extended critical phases 
there is no relevant field in the same symmetry sector as the 
ground state, which is a requirement for the phases to be sta- 
ble. This notion of topological stability will be explained in 



more detail in the section VI dealing with the anyonic spin- 



1/2 chains, where we show in detail that the critical behavior 
of those chains is protected by the topological symmetry. 

As was explained above, there is a relevant operator with 
the same quantum numbers as the ground state at the super- 
conformal point. This operator drives the transition from the 
superconformal point to the Haldane gapped phase on one side 
of the phase diagram, and the extended critical region with Z 2 
sublattice symmetry on the other side. 
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gapless theory 


coset description 


k 


central c 

fc = 5 


large 

k = 7 


SU(2) (k -> 00) 


Zk phase 
Z2 phase 
Zz phase 

superconformal point 


su(2) fe /u(l) 2 fc 
su(2) fe _i x su(2)i/s«(2)& 
su(2) fc _ 4 x su(2) 4 /su(2) k 
su(2) fc _2 x 311(2)2/ 'su(2) k 


c - 2^ 

c — * k+2 

c- 1 6 

(fc+l^fc+2) 

C ~ 2 (fc-2)(fe+2) 
„ _ 3 12 
— 2 fc(fc+2) 


c = 8/7 
c = 6/7 
c = 6/7 
c= 81/70 


c = 4/3 
c = 11/12 
c = 22/15 
c = 55/42 


c = 2 
c = 1 
c = 2 
c = 3/2 



TABLE I: Critical theories in the su(2)fc spin-1 chains for k > 5. 



E. The gapped Haldane phase 



2. Phase boundaries 



In addition to the gapless phases that were discussed in de- 
tail in the previous section, the spin-1 anyonic chains also ex- 
hibits a gapped phase, as can be seen in Figure [4] The proper- 
ties of this gapped phase are strikingly similar to the proper- 
ties of the Haldane phase in the ordinary bilinear-biquadratic 
spin-1 chain. For instance, the point #2,1 = allows for a 
straightforward generalization of the AKLT point of the ordi- 
nary SU(2) model. At this AKLT point, the degenerate ground 
states can be constructed explicitly (see section HIE 3 1. In 



section |HIE4| we discuss the ground states of the open chain, 
and find the degeneracy of the anyonic spin-1 chain can be 
understood in a similar way as the degeneracy of the SU(2) 
model at the AKLT point. Before we deal with the ground 
states at the AKLT point, we first discuss the energy spectrum 
and the phase boundaries of the Haldane phase. 



1. Energy spectrum 

The energy spectrum in the gapped phase is shown in Fig- 
ures [5] and [6] for coupling parameter 9 = 0. It can be seen 
that there exists a quasiparticle band whose qualitative shape 
is identical to the magnon band of triplet excitations of the 
ordinary AKLT point. The complete spectrum is shown at an- 
gle 2 ,i = 0: the ground states occur at momentum K = 0, 
and there exists a quasiparticle band (shown in blue color) 
and a continuum of scattering states (shown in gray color). 
The quasiparticle band is also displayed for coupling param- 
eters #2,1 close to $2,1 = (in red for # 2 .i > 0, in green 
for 6*2,1 < 0). It can be seen that when approaching the crit- 
ical phase with Z3 sublattice symmetry - i.e., for increasing 
9 > - the minimum of the quasiparticle band moves away 
from K — tt towards K = 2ir/3 and K = 4ir/3. When 
decreasing the angle #2,1 < 0, the quasiparticle band remains 
at momentum K = tt, which is consistent with the Z2 sub- 
lattice symmetry of the superconformal critical point. From a 
finite-size scaling analysis of the energy spectra, we confirm 
that the gapped phase does indeed extend over a finite range 
of coupling parameters 9. Figs. [5] and [6] show that the size of 
energy gap (at 6*2.1 = 0) increases from AE(k = 5) f=s 0.16 
to A£(fc = 7) w'0.24. 

This behavior suggests that the qualitative shape of the en- 
ergy spectra at the AKLT point is preserved for all k with the 
energy gap at 6*2,1 = approaching AE(k —> 00) « 0.41^1 



The Haldane phase and the su(2)fe_! x su(2) 1 / su(2) ^ crit- 
ical phase are separated by a superconformal critical point, 
which is located at coupling parameter 6*2,1 s» — 0.1 97r for 
both k = 5 and k = 7. This is very close to the position of the 
phase transition where the Haldane gapped phase gives way 
for a different phase in the ordinary SU(2) spin-1 chain (see 
the phase diagram in figure|4ji, namely 6*2,1 = — arctan(2/3). 

The position of the phase boundary at the other end of the 
gapped phase clearly depends on the level k. Comparing the 
position of this point for k = 5 and k — 7 suggests that it 
moves towards 6*2,1 = i"/2 for increasing k. This scenario is 
consistent with the ordinary model, as can be seen by com- 
paring the phase diagrams of the anyonic and ordinary SU(2) 
spin-1 chain (Figs.|4]and|3]respectively). 



3. Ground states in the periodic chain 
(anyonic equivalent of AKLT point) 

In the ordinary SU(2) spin-1 chain, there exists a point 
within the Haldane gapped phase - the so-called AKLT^point 
- where the ground state can be obtained exactly. At the AKLT 
point, the Hamiltonian penalizes two neighboring spins who 
are in the spin-2 channel. To construct the ground state, it is 
helpful to think of the spin-l's as composed of two spin-l/2's 
which are projected onto the spin-1 channel. In the ground 
state, each of these spin-1/2 forms a singlet with a spin-1/2 
particle that is associated with a neighboring spin-1, as de- 
picted in Figure [7] In this situation two neighboring spin-l's 
will never combine into an overall spin-2 and, therefore, the 
state has zero energy. It can be shown that for periodic bound- 
ary conditions this ground state is non-degenerate^. 

At the corresponding point (angle 6*2,1 = 0) in the phase di- 
agram of the anyonic chains, the Hamiltonian (Eq.( 9|12 1) pe- 
nalizes two neighboring anyons to fuse in the spin-2 channel. 
As for the ordinary SU(2) quantum spin model, the ground 
state can be obtained exactly at this point. In contrast to the 
SU(2) case, there exists a topological symmetry which dic- 
tates that the ground state is degenerate even in the case of pe- 
riodic boundary conditions (we will deal with the open chain 
in the next subsection). One of these degenerate ground states 
is easily found, while the others can be obtained by making 



use of the topological symmetry operator (see section II A for 
details). 

We will present the simplest case of k = 5 here, and give 
the results for arbitrary k in appendixfD] We start by construct- 
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ing one zero energy ground state. For k = 5, the allowed spins 
are 0, 1, 2, and the fusion rules read 



0x0 = 



0x1 
1 x 1 



0x2 = 2 

1x2=1+2 

2x2=0+1 



In particular, the fusion rule 2x1 = 1 + 2 im- 
plies that in the labeling of the Hilbert space, the assign- 
ment (xi-i, Xi, Xi + i) = (2,2,2) is allowed. In addition, 
(aj,_i, Xi, Xi+i) = (2, 1, 2) is allowed as well. Fixing Xi_i = 
Xi + i = 2, one finds that the allowed values of Xi in the trans- 
formed basis are = 0, 1, because and 1 are the two possi- 
ble fusion outcomes of 2x2 = 0+1. Because at the AKLT 
point, only the value xi = 2 is penalized, it follows that the 
state \v ) — |2, 2, . . . , 2) is a zero energy ground state (recall 
that that Hamiltonian is a positive sum of projectors). 

By employing the topological symmetry operators Y;, with 
I = 1,2, we can construct other zero energy ground states. 
The operators Y/ commute with the Hamiltonian, thus the 
states \vi) = Y\ |«o) and \v 2 ) = Y 2 |«o) a ls° nave zero en- 
ergy. It turns out that \ vq) is neither an eigenstate of Y\ nor of 
Y 2 . As a result, the number of ground states is three, which is 
in accordance with the number of particle types in the model. 
We note that Yq is the identity operator. 

The explicit form of the states \vi) and j^) is easily written 
down. First of all, the only basis states with non-zero coeffi- 
cient in \vi) have Xi = 1,2, for all i. Similarly, the only 
basis states with non-zero coefficient in \v 2 ) have Xi — 0, 1, 
for all i. To specify the coefficients, we introduce the nota- 
tion #Z which denotes the number of i's such that xi = I. In 
addition, #(Z, m) denotes the number of i's such that Xi — I 
and Xi+i = in, where we use periodic boundary conditions, 

X L = X . 

Then, we have 

K) = ^2 fi{{x i \)\x ,x 1 ,...x L _ 1 ) (23) 

ir,e{l,2} 

„„ r T 10 <?#£== #(2,l) + #(2,2) 

h{{ Xl }) = {-i)* 2 d- L 4 /2 d\-^d 2 2 

as well as 

|«a) = X] f2({x l })\x ,x 1 , . . .XL^i) (24) 

f2({Xi}) = (-l)#^- i/2 4 /2 #^ #1 • 

Here, d\ and d2 are the quantum dimensions of particles 
with spin-1 and 2 respectively, and are given by d\ = 1 + 
2sin(37r/14) and d 2 — 2cos(7r/7), respectively. 

We labelled the ground states at the AKLT point by \v\) 
with I — 0, 1, 2 for a good reason. In section II A we ex- 



plained that the topological symmetry operators Y\ effectively 
'add' or fuse a particle of type I to the fusion chain. At the 
AKLT point, this notion becomes very explicit. The states \v{) 
are thought of as states of the chain in the I sector. In partic- 
ular, \vq) corresponds to the identity sector. Adding a particle 
of type I, i.e., acting with the operator Yi, gives rise to a state in 



FIG. 7: The AKLT construction of the valence-bond-solid state on 
a finite chain of spin-1 degrees of freedom. Each filled circle repre- 
sents a spin-1/2 variable, each dotted ellipse corresponds to a spin-1 
particle, and and each line connecting two spin-1/2 variables sym- 
bolizes a singlet bond. 



sector I, or \vi) = Yi \vq). Moreover, if one acts with Y\ on the 
state \vj), one obtains a combination of states, which is given 
by the fusion rules. In particular, Y\ \v\) = \vq) + \vi) + \v2), 
Y x \v 2 ) = \vi) + \v 2 ) and Y 2 \v 2 ) = \v Q ) + \vi). Thus, 
loosely speaking, the ground states of the periodic anyonic 
spin-1 chain at the AKLT point form a representation of the 
fusion algebra su(2)fc. Because the modular S matrix diag- 
onalizes the fusion rules, one can easily write down combi- 
nations of the ground states which are also eigenstates of the 

operators Y u namely I^aklt,;) = J2j=o S i-J \ v j)< where Si,j 
is the modular S matrix for (the integer sector of) su(2)s, and 
the sum is over integer values. 

For the explicit form of the AKLT ground states in the gen- 
eral case su(2)fe, we refer to appendix [D] 



4. Ground states in the open chain 
(anyonic equivalent of AKLT point) 

Before describing the structure of the ground states of the 
open anyonic chains at the AKLT point, we briefly review the 
physics of the valence bond solid ground state at the AKLT 
point (6*2,1 = in phase diagram Fig. |4ji of the ordinary 
bilinear-biquadratic spin-1 Heisenberg chain 1 21 . The Hamil- 
tonian at 6 = consists only of the projector onto a total spin- 
2 of two nearest-neighbor spins with a positive sign. Thus, in 
the ground state, a total spin-2 of two-nearest-neighbor spins 
is suppressed. In the usual tensor product basis of local (site) 
states, the valence bond solid ground state is given by 



\^ab) = e 



_ c b 1 a 2t -b 2 a 3 b L ^xa L 



. £ "^-^ \lpab 1 }® \lf>a 2 b 2 ) ® .»® \i> at .b) , 

(25) 

where the summation over repeated upper and lower indices 
is assumed. The local spin-1 state |V>a&) is represented as the 
symmetric part of the tensor product of two spin-1/2 vari- 
ables: 



\lpab 



1 



V2 



{\lpa) ® \lpb) + ® IVO) > (26) 



where ipa denotes one of the two eigenstates of the S z spin- 
1/2 operator, which we label by a — 1,2. The antisymmetric 
tensor e ab enforces a singlet bond of the spin-1/2 variables 
ai + i and 6;. Therefore, the total spin of the two nearest- 
neighbor spin-1 variables, consisting of four spin-1/2 vari- 
ables which are labeled by a;, 6;, a;+i, can only assume 
the values or 1. For a chain with open boundary conditions 
(see Fig. |7|i the first and the last spin-1/2 variables indexed 
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by a-y and do not form a singlet bond. These two spin- 
1/2 variables can add up to a total spin or a total spin 1, 
giving rise to a four-fold degeneracy for the spin-1 bilinear- 
biquadratic chain at the AKLT point with open boundary con- 
ditions. 

With the results above in mind, we will now consider the fu- 
sion basis of the anyonic spin-1 chain, as shown in Figure [T] 
We consider a chain of length L with open boundary condi- 
tions in the sense that variables xq and x L+ i form the ends 
of the chain. In analogy with the above discussion, we as- 
sume that variables xq and xl+i can add up to a total spin of 
xq x xl + i of or 1 in the zero-energy ground states. 

For a given choice of xq and iel+i, we expect that there 
are no zero-energy ground states if \x L+ i — Xq\ > 1 because 
the fusion product xo x x L+ i does not contain nor 1 in this 
case. We expect one ground state to be present if xq x xl+i 
contains or 1, but not both. Finally, if both and 1 appear 
in the fusion product x x x L+ i, we expect two zero-energy 
ground states. There is no S z quantum number in anyonic 
spin chains associated with the 'spins', and the state with total 
spin-1 (or better, topological charge 1) is thus not degenerate. 

The analysis of the previous subsection is helpful in un- 
derstanding the above discussed results. We found that the 
ground states of the periodic chain have a particular form; 
namely, the only basis-states which have non-zero coefficients 
in these states are such that all the Xi take at most two values 
that have to differ by one. Thus, there is a ground state with 
all the Xi £ {0, 1}, one ground state with the Xi £ {1, 2}, etc. 
In addition, the state with all Xi = (k — l)/2 is also a zero 
energy ground state. 

The ground states of the open chain must be such that 
the bulk part of these states does not give an energy con- 
tribution. Thus, for a particular choice of boundary condi- 
tions xo and xl+i, one can construct one ground state if 
\%o — x l+i \ = 1, because there is exactly one corresponding 
zero energy ground state with periodic boundary conditions. 
For Xo — xl + i — 0, there is also one zero energy ground 
state, while for xq = xl+i > 0, there are two zero energy 
ground states. For \xq — x L+ i\ > 1, one finds that there are 
no zero energy ground states. All of this is in accordance with 
the considerations above. 

We computed the ground state degeneracies for all possi- 
ble choices of fixed boundary occupations (xo, xl+i) for both 
the k — 5 and the k = 7 model, and find that the above de- 
scribed picture is indeed the appropriate one. At the AKLT 
point #2,i = 0, the ground state energy is independent of 
the system size. In the Haldane gapped phase away from the 
AKLT point, the ground state degeneracy is not exact and fi- 
nite size effects occur. In Fig. [8] we show the lowest energies 
AEi(x ,x L+ i) = Ei(x Ql x L+ i) - E (x ,x L+ i), i > 1, of 
the su(2)5 spin-1 chain at coupling parameter 6*2.1 = — O.QI71". 
The energy Eq(xq, Xl+i) is the lowest energy of the open 
chain with fixed boundary occupations x and x L+ \, and it is 
not necessarily a ground state energy. By this we mean that 
the state is not a perturbation of a zero energy ground state at 
02 i = 0. For the boundary condition xq = 0, x L+ \ = 2, 
the lowest energy Eq(0, 2) is not a ground state (in the above 
sense) since both A£i(0, 2) and Ai?2(0,2) approach zero 
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FIG. 8: The eigenenergies AEi(xo, xl+i) ■= Ei(xo,XL+i) — 
Eo(xo, xl+i) (i > 1) of the su(2)5 anyonic spin-1 chain with fixed 
boundaries xq and xl+\ as a function of 1/L at 62.1 = — 0.017T. 
The legend at the lower left side indicates the values of xo and 
xl+i- The energy Eo(xo, xl+i) is the lowest energy and not nec- 
essarily a 'ground state energy'. For xo — xl+i = 1, and for 
xq — xl+i = 2, there are two almost degenerate zero-energy states, 
and AEi(xq,xl+i) corresponds to the finite-size splitting of the 
two ground states that decays exponentially with system size (see 
the inset). 



in the limit 1/L — > 0, as demonstrated in Fig. [8] For the 
boundary condition xo = 1, xl+i = 1, as well as xq = 2, 
xl + i = 2, the ground state is two-fold degenerate, and the 
splitting of the two ground state energies at finite system size 
L decays exponentially in 1/L, as illustrated in the inset of 
Fig. [8] Again, this is in agreement with the above discus- 
sion because lxl = 0+1 + 2 and 2x2 = + 1 (for 
su(2)s), i.e. both fusion products allow for a total spin and 
a total spin 1, For all remaining possible boundary condi- 
tions, there is one ground state, as can be seen from Fig. [8] 
where AE^Xq, x l+1 ) approaches a finite energy in the limit 
1/L — > 0. We also verified this scheme for the su(2)7 model, 
and for different values of 6*2,1 in the gapped phase. 



IV. ANYONIC SU(2) fe SPIN-1 CHAINS: EVEN k > 6 

In the previous section, we discussed in detail the odd-k 
anyonic spin-1 chains. We found that the phase diagram of 
these models (see Pig.[3]», bears great resemblance to the phase 
diagram of the SU(2) spin-1 chain (see Fig.Q. We observed 
one striking difference between the ordinary and the anyonic 
spin-1 chains; namely, the absence of a (gapped) 'dimerized' 
phase in the case of the anyonic spin-1 chains. In this section, 
we present our result for the even-fc anyonic spin-1 chains. 
For even k, the phase diagram is very similar to the case of 
odd k with the exception of an additional gapped phase which 
resembles the dimerized phase of the SU(2) spin-1 chain. 

In this section, we focus on the case k = 6; however, our 
analysis for k — 8 indicates that the case k — 6 is generic 
for k even. The generic structure of the phase diagram for 
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su(2)k spin- 1 chains (even k) 




FIG. 9: (color online) Phase diagram of the even-fc anyonic 
su(2)fc spin-1 chain in a projector representation |9]( where Ji = 
— sin(#2,i) and J2 = cos(6*2,i). The locations of the phase bound- 
aries correspond to the case k — 6. Some of the phase boundaries 
move with increasing (even) k; the arrows indicate the direction of 
the change. 



even fc > 6 is analogous to the generic structure of the phase 
diagram for odd k > 5. We note that the case k = 4 is special 
and will be considered in detail in the following section. 

The fact that the phase diagrams for k even and odd dif- 
fer is a very interesting feature of our model. As far as we 
are aware, this is the first time that a dependence on the par- 
ity of the level k has been observed. As we will point out in 
the discussion, Koo and Salem 40 considered a closely related 
loop model which contains a continuous parameter that plays 
the role of the discrete level k. The model considered by Koo 
and Saleur does not show any sign of the 'even-odd' effect we 
observe. It would be very interesting to understand the differ- 
ences and similarities of these two models in greater detail. 

The phase diagram of the k = 6 anyonic spin-1 chain is 
presented in Figure [9] We will discuss the similarities and 
differences of this phase diagram to the phase diagram of the 
case k = 5 (Fig. pjY The locations of the phase boundaries in 
Figure [9] correspond to the case k = 6. As was the case for 
k odd, we observe that some of the phase boundaries change 
upon increasing the value of (even) k. The direction of the 
movement of the phase boundaries is indicated by the arrows 
in the phase diagram. 

Comparing the phase diagrams for odd and even k in Fig- 
ures|4]and|9j we first note that large parts of the phase diagram 
have a similar structure. At angle 6* 2 ,i — 0, we encounter a 
gapped Haldane phase, precisely as in the case of odd k. At 
angle 62,1 ~ — 0.197T, there is a phase transition that is de- 
scribed by a N = 1 supersymmetric minimal model from the 
Haldane phase into an extended critical region (we will com- 
ment on the latter critical region below). At the other end 
of the gapped Haldane phase, there is a phase transition at 
angle #2.1 ~ 0.097T (for k = 6) to a critical region that ex- 
hibits a Z 3 sub-lattice symmetry and is described by the coset 
su(2)4 x su(2)k-i/ 'su(2)k (we note that the corresponding 



critical region for odd k is described by the same CFT). This 
critical region extends all the way to 6*2. 1 = 37r/4 at which 
point there is a first order transition to a critical region with 
Z% sublattice symmetry. So far, the phase diagram for even k 
has the same structure and phases as the one for odd k. 

The phase diagrams for odd k versus even k begin to di- 
verge at the angle where for k odd, the critical region with 
Zk sublattice symmetry transitions to a critical phase with Z2 
sublattice symmetry. While the former (Zf.) critical phase 
also appears for k even, the latter (Z2 critical phase) does 
not; rather, there is a phase transition at 6*2.1 ~ 1.417T (for 
k = 6) to a gapped phase. This gapped phase is char- 
acterized by broken translational invariance, as signified by 
a zero-energy ground state at K = ir present at the angle 
62,1 = 37r/2. In addition, there are (fc + 2)/2 degenerate 
ground states at momentum K = with topological quantum 
numbers (0, 1,2, ... , fc/2). The zero energy ground state at 
K = it is in topological symmetry sector fc/4. Clearly, the 
nature of this 'dimerized' gapped phase differs from the Hal- 
dane gapped phase. 

Between the 'dimerized' gapped phase and the Haldane 
gapped phase, we find an extended critical region. Due to 
the rather small extend of this critical region and the fact that 
we could not study large enough systems (the dimension of 
the Hilbert space increases with k), we have not been able to 
determine which CFT describes this extended critical region. 

It is interesting to note that the structure of the phase dia- 
gram for even k bears closer resemblance to the phase diagram 
of the SU(2) bilinear-biquadratic spin-1 chain, (see Figure [3]) 
than to the phase diagram for odd-fc anyonic spin-1 chains. In 
particular, both the phase diagrams of the ordinary SU(2) spin- 
1 chain and the even-fc anyonic spin-1 chain exhibit dimerized 
phases in the area surrounding the angle 6*2,1 = 37r/2. It ap- 
pears that for increasing even fc, the phase diagram of the any- 
onic chain gravitates towards the phase diagram of the SU(2) 
chain. Our results for the fc = 8 anyonic chain are consistent 
with this picture. 

The phase diagram for the fc = 6 anyonic spin-1 chain dis- 
plays a unique feature; namely, its structure is symmetric in 
the line through the points 6*2,1 = 37r/4, 77r/4. The under- 
lying reason is that the fusion rules of the su(2)j. theory are 
symmetric under the exchange j fc/2 — j, where the la- 
bels j take the values j = 0, 1/2, ... , fc/2. In the case of 
fc = 6, this symmetry exchanges anyon spins 1 <-> 2. The 
location of the symmetry points follow from our parametriza- 
tion of the hamiltonian, as given in equation |9|. We point out 
that this symmetry only relates the sets of energy eigenvalues, 
but not the possible degeneracy of the levels or their angular 
momenta. 

For example, the energy levels levels at the point 6*2.1 = tt - 
where the system is described by the Z 6 parafermion theory - 
are identical to those at angle 6*2.1 = n/2. At the latter point, 
the system is described by the coset su(2)2 x su(2)4/su(2)g, 
which for fc = 6 corresponds to the Zq parafermions. We note 
that the momenta of the states are not identical. 

Similarly, the energies of the levels in the dimerized gapped 
phase are the same as the energies of the levels in the Haldane 
phase, even though the nature of these gapped phases is very 
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different. We will return to this issue below. Finally, we note 
that the phase transition from the dimerized phase to the crit- 
ical region in between the dimerized phase and the Haldane 
phase is given by an N = 1 supersymmetric model. As far 
as we can tell from our numerics, this is only true for the case 
k = 6. For k = 8 and higher, we have not been able to deter- 
mine the CFT describing this phase transition. 



V. ANYONIC SU(2) fc SPIN-1 CHAINS: k = 4 

Having discussed the anyonic spin-1 models for odd k > 5 
and even k > 6, we finally turn our attention to the remaining 
case k = 4. We already pointed out in the introduction that 
the phase diagram for k = 4 has a different structure than the 
phase diagrams for other values of k. The underlying reason is 
that the spin-1 particle is special in this case. The symmetry of 
the fusion rules under the exchange j O k/2 — j implies that 
j = 1 is mapped onto itself for k — 4. In addition, k = 4 is 
the lowest k for which a general fusion rule 1x1 = 0+1 + 2 
applies. We refer to the discussion in section VII for more 
details. 



A. Hilbert space and Hamiltonian 

The basis of the su(2)4 spin-1 chain is depicted in Fig. [T] 
Each labeling {xj}j = o,....z,-i € {0, ^, 1, §,2} that satisfies 
the fusion rules at the vertices corresponds to a different basis 
state. In fact, the Hilbert space of the su(2)4 spin-1 chain splits 
into two independent sectors: the fusion rules impose that the 
local basis elements are either all integer valued or all half- 
integer valued. We shall use the following terminology: 



Integer sector (IS): {xi}i 



€ {0,1,2}. 



Half-integer sector (HIS): {xi}i 



■L-i e ih §} 



We will only consider periodic boundary conditions for the 
su(2)4 chain, i.e., xl = xo. 

We find that the differences in behavior between the IS and 
HIS su(2)4 spin-1 chains are rather subtle. We will first de- 
scribe the behavior of the model in the IS sector, followed by 
a discussion of the HIS sector. 

As a first minor difference, we note that the number of states 
in the HIS is given by 2 L + 5l,o, where L is the length of the 
chain. In the IS sector, however, the number of states is 2 L + 1 
when L > is even and 2 L — 1 when L is odd. The additional 
state in the even-L IS occurs at momentum K = vr, while 
the additional state in the odd-L HIS occurs at momentum 
K = 0. Those are the only differences; the remaining 2 L 
(2 L — 1) states where L even (odd) have the same momenta in 
the integer and half-integer sectors. 

As we did for k > 5, we represent the Hamiltonian of the 
su(2)4 spin-1 chain in terms of the projectors onto the 1 and 2 
channels with couplings Ji and J 2 , respectively. These cou- 
plings are parametrized by an angle 02, l where J 2 = 0036*2.1 



and Ji = — sin fe, 1. Explicitly, the Hamiltonians read 

H is =4) = E cos 9 ^ P i% - sin °^ P < % ( 27 > 

i 

#his 4) = E cos e *,i P 8ma 2) ii^ IS • (28) 



The explicit form of the projectors are given in ap- 



pendix C 3 a 



B. Phase 



diagram in the integer Hilbert space sector (IS) 



The phase diagram of the IS su(2)4 spin-1 chain (Hamil- 
tonian given in eq. (|27]i) is shown in the left most panel of 



Fig. 10 The phase diagram consists of two extended gapped 
phases which are separated by two extended gapless regions. 
The two phase transitions between the gapped phase with a 
Z 3 -sublattice structure and the two gapless regions are first 
order. However, the phase transitions into the gapped phase 
with a ^2-sublattice structure are continuous. 

The critical behavior of the gapless regions is described by 
the Z 2 orbifold theory of the u(l)-compactified boson with 
central charge c = 1. Interestingly, the compactification ra- 
dius varies continuously as a function of 02, 1 in the gapless 
regions. We found it difficult to determine the range of com- 
pactification radii which are realized in the model. The reason 
is that the finite size data makes it difficult to determine the 
location of the transition between the gapped phase with the 
Z 2 -sublattice structure and the critical regions. We will de- 



vote a separate subsection V D to the issue of the location of 
these phase boundaries, dealing with the IS and the HIS at the 
same time. 



1 . Gapped phases (IS) 

Gapped phase 2j i = it/2 - The gapped phase containing 
the point 02. 1 = tt/2 extends from 02. 1 = to 02. 1 = 37r/4. 
These phase boundaries are easy to locate because the tran- 
sitions to the gapless regions are first order transitions, as we 
will show in section [VDl 

This gapped phase has a Z3 sublattice symmetry, which re- 
sults in a three-fold degenerate ground state for system sizes 
that are a multiple of 3. These ground states occur at momenta 
K = 0, 27r/3,47r/3 and their exact form can be established 
throughout the whole gapped phase. 

At angle 02 1 = tt/2, the Hamiltonian can be solved ex- 
actly. At this point, the Hamiltonian reduces to the equal sum 
of two projectors, namely onto the spin-0 and spin-2 channels 
(in addition, there is also a constant term — L). Throughout 
the region < 02,i < 37r/4, the Hamiltonian is a sum of two 
projectors with positive coefficients (for a matrix representa- 



tion of the Hamiltonian see appendix C 3 1. The three degen- 
erate ground states are build from the basis state of the form 
I llai lla2 Ila3 * * * ) H fl L/3) an d its two translations, where 
the a-i represent the states \a) 3i = (|0) — \2))/^/2 at position 
3i. These three states can easily be combined to form three 
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su(2)4 spin-1 chains 

a) integer sector J1-J2 b) integer sector J0-J2 c) half-integer sector J1-J2 




O Kosterlitz-Thouless O Ising 2 O Parafermions • 4-state Potts O super CFT 

FIG. 10: (color online) Phase diagrams of the su(2)4 spin-1 chain in the integer sector and half-integer sector and different projector represen- 
tations. The colored circles indicate special points in the c = 1 gapless phase that can be matched to the labelled conformal field theories. 



momentum eigenstates. These ground states have energy — L 
and are eigenstates of the two projectors with eigenvalue 0. 
The latter explains that these ground states persist throughout 
the whole gapped phase. 

Gapped phase 6*2,1 = 37r/2 - In the gapped phase sur- 
rounding the point 62,1 = 3tt/2, the spectrum exhibits a 
Z2 sublattice symmetry and a cosine-shaped quasi-particle 
dispersion. For even-L system sizes, the ground state 
is threefold degenerate, with the ground states occur- 
ring at momenta K — 0,0, ir. Two of the three ground 
states at angle 6*2.1 = 37r/2 consist of states of the form 
I62 I63 ■ ■ ■ , I&L/2) an d hs translation by one site, where 
|fo) 2 i = (|0) + \2))/y/2 at site 2i. These two state can 
be combined to form the two ground states at momenta 
K = 0, 7t. The state | 111 ... 1) is the third ground state and 
has momentum K = 0. For odd system sizes, this state is the 
only ground state of the system. 



2. Gapless phases (IS) 

The critical behavior of the su(2)4 spin-1 chain is particu- 
larly interesting. We find that the critical behavior depends 
continuously on the angle describing the interaction. At par- 
ticular values of the angle O2, 1, the behavior matches particu- 
lar CFTs with central charge c = 1. In particular, these CFTs 
are the Z2 orbifolds of a boson compactified on a circle of 
radius R = \/2p : For p integer, these are rational co nform al 
field theoriegSES, described in detail in the appendix 



E5 



In 

this section, we will limit the discussion to the most prominent 
features of these theories. In section VD we will point out 
the particular orbifold theories that are realized in the su(2) 4 
spin-1 chain. 



To identify the critical theories describing the critical be- 
havior as a function of the angle, we employ the standard 
technique of first shifting the spectrum such that the ground 
state has zero energy, followed by a rescaling of the energy to 
elucidate the conformal nature of the spectrum. 

By means of this procedure, we identified several of the 
c = 1 orbifold theories. These theories 'share' several opera- 
tors that appear in the spectrum throughout the critical region. 
These operators are the ground state with h = h = 0, two 
twist fields (T12 with scaling dimension h a + h a = 1/8, two 
twist fields n,2 with scaling dimension h T + h T = 9/8, a field 
with scaling dimension 2, and finally, two fields $12 with 
scaling dimension p/2. For p = 1, the fields just described 
exhaust the full list, but in general, there are p — 1 additional 
fields <p\ with scaling dimension These fields, as well 
as the associated momenta and topological symmetry sectors 
are given in table [TT] We checked that the assignments of the 
topological symmetry sectors are compatible with the fusion 



rules of the orbifold CFTs (for details, see appendix E 5 1 



For various values of p, the orbifold theories are also known 
under specific names, such as the Kosterlitz-Thouless theory 
(p = 1), the theory of two decoupled Ising models (p = 2), the 
Z4 parafermion CFT (p = 3), the 4-state Potts model (p = 4), 
and the superconformal minimal model with c = 1 (p — 6). 

We identified several of the c = 1 Z2 orbifold theories, 
including the ones with p = 1,2,3,4,5,6. In the left side 



of Figure 1 1 we show the energy spectra associated with the 
p = 1 and p = 3 orbifold theories in the Z2 critical region. 



In the right side of Figure 11 we display the energy spectra 
associated with the p = 2 and p = 6 orbifold theories in the 
Z4 critical region. 

In table [Til] we list the locations of some of the critical 
points as extracted from the numerical data. The procedure we 
followed to obtain these locations will be described in more 
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TABLE II: The scaling dimensions (h + h) of the operators of the Z2 orbifold of the compactified boson on a circle of radius R = ^J2p for 
some integer p. The following abbreviations are used: sCFT = the (minimal) superconformal CFT with central charge c = 1, Potts = 4-state 
Potts CFT, pCFT = Zi parafermion CFT, (Ising) 2 = square of the Ising CFT, KT = Kosteiiitz-Thouless theory, equivalent to the compactified 
boson theory u(l)g. We also list the numerically observed topological quantum numbers (F-symmetry: 1/0 = 2/2 = 2, yi/2 = 2/3/2 = 0, 
1/1 = — 1) and momentum quantum numbers K at which the fields appear in the various critical regions. The symmetry sectors of the fields 
with scaling dimensions p/2 depend on p. This is a consequence of the fact that the field with scaling dimension (p — l) 2 /2p at radius p 
corresponds to the field with dimension p/2 at radius p — 1. 



p 


theory 


Zi 


z 2 


1 


Kosterlitz-Thouless 


0.755tt* 


-O.OLt 


2 


Ising 2 


0.77tt 


-0.047T 


3 


parafermion 


0.807T 


-0.087T 


4 


4-state Potts 


0.83tt 


-0.137T 


5 




0.88tt 


-0.177T 


6 


superCFT 


0.92tt 


-0.207T 


7 




0.96vr 


-0.23tt* 


8 




0.98vr 


-0.24tt* 



region) is 



TABLE III: The approximate locations of some of the critical theo- 
ries of the su(2)4 spin-1 chain l |27| l in the integer sector (IS) are listed 
for both the Z2 and the Z4 critical regions. The angles without as- 
terisk are obtained directly from exact diagonalization for L — 20, 
i.e. we matched the momentum resolved spectrum to the CFT. The 
angles with an asterisk were obtained by using the relations between 
the angles 62 and 64, as explained in the text. We only list those val- 
ues of p for which we could match the CFT description beyond any 
doubt. 



detail in section |VD| The location 9 ss — 0.207T of case p = 6 
- the superconformal theories - is very close to the location 
of the superconformal point for the su(2)fc spin-1 chains with 
k > 5, namely #2,1 ~ — 0.197T. 

The location of the superconformal point in the Z4 criti- 
cal region is 6*2,1 ~ 0.927T. In general, the relation between 
critical angle in the Z 2 critical region (which we will for now 
denote by 62, similarly, 84 denotes the angle in the Z4 critical 



tan -1 (1 + tan 2 ) 6» 2 = -tan 1 (l + tan0 4 



(29) 



The spectra in Figure 11 illustrate the different sublattice 
symmetry for the two gapless regions. In these spectra, we 
also indicate the topological symmetry sectors of some of the 
low-lying states. In the case of su(2)4, the topological symme- 
try operator Y has three distinct eigenvalues, which are given 
by yo = V2 = 2, 2/1/2 = 2/3/2 = and y x = -1. We will thus 
use the labels y = 0, 1/2, 1 for these sectors. 

The presence of the different critical models with the same 
central charge c = 1 indicates the presence of a marginal op- 
erator that drives the 'transition' between the different critical 
theories and that gives rise to continuously varying critical ex- 
ponents. Indeed, all the orbifold models share a marginal op- 
erator 9 with scaling dimension 2 whose topological symme- 
try coincides with that of the ground state. It is this operator 
which is responsible for the critical region with continuously 
changing exponents. It proved difficult to locate the phase 
transition between the critical regions and the gapped phase 
around 6*2,1 = 37r/2. One reason might be that the transi- 
tion to the gapped phase is also driven by a marginal operator, 
which allows for large finite size effects that thwarts the local- 
ization of these critical points. 
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FIG. 1 1 : (color online) The su(2)4 chain - integer sector: Energy spectra at various in the gapless phases of the phase diagram in Fig.|10[i). 
The energy spectra have been rescaled to match the conformal field theory prediction given in Eq. l |22[ (. Green squares indicate the location 
of the primary fields, red circles the descendant fields. The topological symmetry sector is indicated by the violet index. Data shown are for 
system size L = 20. 



C. Phase diagram in the half-integer Hilbert space sector 
(HIS) 

The behavior of the su(2)4 spin-1 chain in the half-integer 
sector mimics very closely that of the integer sector. The 
phase diagram is presented in the rightmost panel in Fig [10] 
The phase boundaries are located at the same positions, but the 
details of the observed phases differ slightly. In the following 
discussion of the HIS su(2)4 spin-1 chain, we will emphasize 
the differences between the two sectors. 



1. Gapped phases (HIS) 

As was already noted above, there are some differences in 
the dimensions of the Hilbert spaces in the IS and the HIS, 
respectively. As a consequence, the IS and HIS models have 
different sublattice structures in the gapped phases. Namely, 
in the half-integer sector, the ground state occurs at momen- 



tum K = 0, and it is non-degenerate. All other features of 
the gapped phases in the half-integer sector are very similar to 
those observed in the integer sector. 

Gappedphase 6 2 ,i = 7r/2 — In gapped phase that surrounds 
the angle 62,1, the ground state is non-degenerate and occurs 
at momentum K = (there is no sublattice structure). The 
model can be solved at angle #2,1 = tt/2: the ground state 
can be expressed as follows, 

\GS)= (-l) #(3/2 ' 3/2) |V2,...,*L-l>. 

^=1/2,3/2 

Here, #(3/2, 3/2) denotes the number of times the sequence 
(xi, Xi+i) — (3/2, 3/2) occurs in the state \xq, X2, ■ ■ ■ , Xl-i) 
(note that periodic boundary conditions impose xl — Xq). 
As was the case for the ground state(s) at #2.1 = tt/2 in the 
IS, this state is in fact the ground state throughout the whole 
gapped phase, i.e. for angle < 6*2,1 < 37r/4. 

Gapped phase 6 2 .i — 3tt/2 - In the gapped phase that 
surrounds the angle 6*2,1 = 37r/2, the ground state is non- 
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degenerate and occurs at momentum K = (there is no sub- 
lattice symmetry). At Q% \ — 3ir/2, the ground state is given 
by 



\GS) = 



Sj=l/2,3/2 



(-1) 



#(1/2,3/2) 



X ,X 2 , 



■,xl-i) 



All basis states contribute to the ground state. The 
sign of a term is given by the number of times the se- 
quence (xi,Xi + i) = (1/2,3/2) occurs in the basis state 
|xo, X2, ■ ■ ■ , Xl-x) (periodic boundary conditions are as- 
sumed). 



D. The location of the phase boundaries 

To locate the boundaries of the gapped and critical regions 
of the su(2)4 spin-1 chain, we consider the ground state energy 
as a function of the interaction angle 0. The analysis is most 
easily carried out by using an alternative parametrization of 
the Hamiltonian. Two spin-1 anyons can fuse into either a 
spin-0, a spin-1 or a spin-2 anyon; therefore we can write the 
Hamiltonian in terms of projectors onto the spin-2 and spin-0 
channels, instead of the spin-2 and spin-1 channels as we did 
in eq.(([27)). By making use of the relation I = p(°) + pW + 
p( 2 ', we find that the Hamiltonian 



H 



(fc=4) q p (2) . „ p (0) 

J2— JO - C0S - Sill 02,0-^i 



(30) 



2. Gap! ess phases (HIS) 

As in the integer sector, the phase diagram in the half- 
integer sector has two extended regions where the model is 
critical. The criticality is again described by Z 2 orbifold 
models. We identified the orbifold models with parameters 
p = 2, . . . , 9. Some of the critical angles are given in table IV 
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TABLE IV: The approximate locations of some of the critical theo- 
ries of the su(2)4 spin-1 chain in the half integer sector (eq (28}) for 
both the Z2 and the Z4 critical regions. The angles without aster- 
isk are obtained directly from exact diagonalization for L = 20 by 
matching the momentum resolved spectrum to the CFT The angles 
with an asterisk were obtained by using the relations between the an- 
gles 82 and $4, as given in eq. \29\ . We only list those values of p for 
which we were able to match the CFT description beyond any doubt. 



The difference between the two gapless regions in the half- 
integer sector lies in the momentum quantum numbers, as in- 
dicated in Table [II] The topological symmetry sectors in the 
HIS coincide with those found in the IS (see Table |Hj). While 
this is to be expected for topological quantum numbers, it nev- 
ertheless shows that our results are consistent. 

A major distinction between the IS and the HIS phase di- 
agram of the su(2)4 spin-1 chain is the order of the orbifold 
theories. By comparing the leftmost and the rightmost panels 
of Fig. [TO] it can be seen that in the integer sector, the orbifold 
theories appear in ascending p order when moving away from 
the first order transition points, while in the half-integer sec- 
tor, the orbifold theories appear in descending p order when 
moving away from gapped phase I. We exploit this result in 
locating the position of the critical endpoint of one of the the 
gapped phases (see following subsection). 



is related to the Hamiltonian of equation ( |27) 



H 



via 



tan 1 



'2,1 = 



tan 8 2 a 
1 + tan 0. 



(31) 



(32) 



2.0 



up to an unimportant shift in energy. 

The ground state energy as a function of the angle 6*2,0 is 
given in figure 13 for a chain of size L = 18. The kinks in 
the ground state energy indicate that there are two first order 
phase transitions. These first order phase transitions mark the 
boundaries of the gapped phase located at — 7r/2 < 0% q < 
in the new angle variable 02,0 (0 < 6*2,1 < 37r/4 in terms of 
the original variable 02, 1, 



lOi 



see phase diagram in Figure 
To identify the location of the continuous transition be- 
tween the other gapped phase and the neighboring gapless 
phases, we plot the first and second derivatives of the ground 
state energy per site. From these derivatives, it can be con- 
cluded that these transitions are roughly located at 02. = 7r/2 
and 02,o = 7T- In terms of the original variable 02, 1, these 
locations correspond to 02. 1 = — 7r/4 and 02, 1 = 7r. This 
conclusion is corroborated by figure [14] in which we plot the 
ground state energy in the gapped phase surrounding the angle 
02.i = 37r/2, i.e. 02,o = 37r/4 for system sizes ranging from 
L = 8 to L = 20. 

In order to locate the phase boundaries, we also consid- 
ered the structure of the orbifold CFTs describing the gapless 



phases (we refer to Appendix E 5 for more details on the orb 



ifold CFTs). We know that throughout the critical region, two 
fields with scaling dimension h + h = 1/8 and two fields 
with scaling dimension h + h = 9/8 must appear. In addi- 
tion, there are several fields with scaling dimensions A 2 / (2p) 
(A = 1, ... ,p — 1) for some value of p. Depending on the 
sector (IS or HIS), and depending on the critical region, these 
fields appear at different momenta, as detailed in table [II] This 
table also includes our numerical results for the topological 
symmetry sectors of the various fields. 

The structure of the critical theories describing the critical 
region allows us to numerically determine the value of p as 
a function of the angle 6*2,0- Moreover, in doing so, we will 
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FIG. 12: (color online) The su(2)4 chain - half -integer sector: Energy spectra at various points in the gapless phases of the phase diagram 
displayed in Fig.|10|:). The energy spectra have been rescaled to match the conformal field theory prediction given in Eq. {22J. Green squares 
indicate the location of the primary fields, red circles the descendant fields. The topological symmetry sector is indicated by the violet index. 
Data shown are for system size L = 20. 




FIG. 13: (color online) The ground-state energy per site (upper 
panel) and its first and second derivative (middle and lower panel, 
respectively) of the IS su(2)4 chain. Data shown is shown for system 
size L = 18. 



gain insight into the locations of the phase boundaries. We 
proceed as follows: We first shift the spectrum such that the 
ground state has energy zero, and we rescale the spectrum 
such that the two degenerate lowest fields with topological 
eigenvalue y — have energy 1/8. Since these fields are 
always among the low-lying fields, finite size effects are in- 
significant. After shifting and rescaling the energy, we focus 
on the two states corresponding to the fields with scaling di- 
mensions l/(2p) and 4/(2p). By equating the numerical en- 
ergies to the p-dependent predictions from the conformal field 
theory, we obtain a numerical estimate of p as a function of 
the interaction angle. We note one has to be watchful of level 
crossings when using this procedure. 



In Figs. 15 and 16 we display the numerically obtained val- 
ues for p as a function of the angle for system size L = 20. 
In these Figures, we also show the energy of the state corre- 
sponding to the field with scaling dimension 9/8. The range 
of angles 6*2.0 over which the field with scaling dimension 9/8 



is constant is shaded in Fig. 15 the shaded region includes all 



angles for which the energy associated with the field multi- 
plied by eight takes values between 8.9 and 9.1. It is imme- 
diately apparent that the two independent numerical estimates 
of p agree very well in the range < #2,0 < vr/2. This applies 
to both integer and half-integer sector. In addition, the energy 
of the state corresponding to the field with scaling dimension 
9/8 agrees very well with the prediction over this range. Thus 
our numerical data is consistent with the picture that the Z2 
critical region extends over the range < 62, < tt/2, giving 
way at 6*2,0 = 7r/2 to the gapped phase with a Z2 sub-lattice 
structure. 

We do not include similar Figures for the Z 4 critical region, 
but note that they show very similar behavior. This leads to the 
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FIG. 15: Numerical estimate (system size L = 20) of parameter p 
from the eigenenergies that are associated with operators with scaling 
dimensions 1 /2p (red squares) and 4/2p (blue dots) in the IS gapless 
phase with twist operators at k — and k — n. The parameter p is 
about 1 at #2,0 = 0, and grows to about 9 at 82,0 = tt/2. For angles 
#2,0 > 7r/2, the estimates of p obtained from the two operators start 
to deviate. The black dots correspond the scaling dimension of the 
fields with dimensions h = 1/8 and h = 9/8 multiplied by eight, 
as obtained from exact diagonalization. The shaded region indicates 
the range of #2,0 for which the latter dimension lies between 8.9 < 
8h < 9.1. 
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FIG. 16: Numerical estimate (system size L = 20) of parameter p 
from the eigenenergies that are associated with operators with scal- 
ing dimensions l/2p (red squares) and 4/2p (blue dots) in the HIS 
gapless phase. The black dots correspond to the scaling dimension 
of the fields with dimensions h — 1/8 and h = 9/8 multiplied by 
eight, as obtained from exact diagonalization. The shaded region in- 
dicates the range of 82.0 for which the latter dimension lies between 
8.9 <8h< 9.1. 



in the integer and half-integer sectors, respectively, both the 
low- and the high-p orbifold CFTs are observed near the first 
order phase transitions at # 2 ,o = 0. In the IS, we identified 
the p = 1 orbifold CFT (see Figure 111, which suggests that 
the gapless phases in the IS su(2)4 chain include the orbifold 
CFTs starting at the lowest integer value p = 1. In contrast, 
in the HIS, we were able to match the spectrum for the p = 2 
CFT, but we did not find evidence that the p = 1 model exists 
in the phase diagram. Moreover, we found that the values of 
p in the HIS increase quite rapidly when decreasing 62,0 to 
zero. We were able to identify the orbifold CFTS up to p = 9. 
The reason is that the abundance of low-lying (primary) fields 
in the high-p orbifold theories requires large system sizes to 
identify these CFTs with sufficient accuracy. 

Further insight into the critical phases can be gained by 
considering the topological sectors of the various operators 
in the spectra for the integer values of p (see Table [H}. All 
orbifold CFTs (i.e., all p) include a marginal operator with 
conformal dimension h = 2. This marginal operator has mo- 
mentum K = and topological quantum number y , i.e., it 
has the same quantum numbers as the ground state. It is this 
marginal operator which causes the continuously varying crit- 
ical behavior within the gapless phase. With increasing p, the 
number of fields whose scaling dimensions are smaller than 
two increases. However, these fields are not relevant because 
their topological and/or momentum quantum numbers differ 
from those of the ground state. The lowest-p orbifold CFT for 
which there exists an additional marginal operator with the 
same quantum numbers as the ground state is p — 9. For gen- 
eral p, this operator has scaling dimension h + h = 36/ (2p) 
(see Table |H|). The existence of an additional relevant opera- 
tor for p > 9 (as 36/2p surpasses two for p > 9) suggests 
that the range of p values which are realized in our model is 
p = 1, 2, . . . , 9. Moreover, if the marginal operator which 
first appears for p = 9 is indeed the operator which is driving 
the phase transition, it is not surprising that the location of the 
continuous phase transition is hard to determine. 

In conclusion, we provide evidence that the orbifold CFTs 
with p = 2, 3, . . . , 9 are realized in both IS and HIS of the 
su(2)4 spin-1 anyon chain, while the p = 1 orbifold CFT ap- 
pears only in the IS. We note that the su(2)4 anyonic spin-1 
chain has some similarities with the one-dimensional quan- 
tum Ashkin-Teller model 43 . The one-dimensional quantum 
Ashkin-Teller model, which is an anisotropic version of the 
two-dimensional Ashkin-Teller modeF^ also has a line of 
critical points on its self-dual line, realizing the orbifold CFTs 
with p = 1,2, 3, 4, in addition to two gapped phases, one of 
which has a Z2 sub-lattice structure. 



location of the boundaries of this critical region being #2,0 = 
7r and #2,0 = 37r/2, the latter being the location of the first 
order transition. 

We also studied the various values of p which are realized 
in the su(2)4 spin-1 model. The region close to the first order 
transitions is most suitable for identifying the various orbifold 
models because of drastic changes in the spectrum in this re- 
gion. Since the Z 2 orbifold CFTs appear in opposite order 



VI. ANYONIC SU(2) fc SPIN-1/2 CHAINS 

In this section, we discuss the results of our study of the 
su(2)fe spin-1/2 anyonic spin chains for k = 2,4, 5. The case 
k = 3 is the original 'Golden Chain' model, which marked the 
beginning of the study of anyonic quantum spin chains^. In 
the latter publication, it was established numerically that for 
both antiferromagnetic as well as ferromagnetic interactions, 
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FIG. 17: (color online) The su(2)2 spin- 1/2 chain: Energy spectra have been rescaled to match the conformal field theory prediction given in 
Eq. \22) . Green squares indicate the location of the primary fields, red circles specify the descendant fields. The topological symmetry sector 
is indicated by the violet index. Data shown is for system size L = 40. 
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the system is critical, and that the system can be described by 
the tri-critical Ising model and the Z3 parafermion eft (three- 
state Potts model criticality), respectively^. 

In addition, it was shown that the model can be mapped 
onto an exactly solvable model, namely a particular 'restricted 
solid-on-solid' (RSOS) model 16 . This mapping is applicable 
to arbitrary k, and thus the critical behavior of the spin-1/2 
anyonic chains is described by the fc-critical Ising model for 
AFM interactions and Z^-parafermions for FM interactions 16 . 

Finally, it was conjectured in Ref. ([I6i l that the critical- 
ity of these spin-1/2 anyonic chains is not merely due to a 
fine tuning of parameters, but is in fact protected by a non- 
local, topological symmetry of the model. This implies that 
the model remains gapless if a perturbation which preserves 
both the spatial and the topological symmetry is added to the 



model. This property is essential for the nucleation of a new 
topological liquid as a result of interactions between anyons^. 

In this section, we consider the topological symmetry prop- 
erties of the su(2)j. spin-1/2 chains and explain why the crit- 
icality is topologically protected for all finite k. Explicit 



Hamiltonians are given in appendix C 2 



The numerically obtained spectra for both AFM and FM in- 
teractions are given in Figures 17 18 and 19 for k = 2, k = 4 
and k = 5. The spectra were obtained by exact diagonal- 
ization of the Hamiltonian, followed by shifting and rescaling 
of the spectrum in order to match the conformal field theory 
predictions. 

The numerical results confirm that the spin-1/2 su(2)fc 
chains are described by the fc-critical Ising model for AFM 
interactions and the Zk -parafermion CFT for FM interactions. 
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Details of these CFTs are given in appendices E 1 and E 4 

In the remainder of this section, we discuss the assignment 
of topological symmetry sectors to the states in the energy 
spectra, as indicated in Figures 17 18 and 19 The topologi- 



cal symmetry sectors were obtained by acting with the opera- 
tor Y on the eigenstates. Because Y commutes with both the 
Hamiltonian and the momentum operator, and because most 
states are non-degenerate, it follows that the eigenstates of the 
Hamiltonian (in the momentum representation) are also eigen- 
states of the topological operator Y, 

We begin the analysis with a general observation. A topo- 
logical symmetry sector is assigned to each state in the spec- 
trum. Moreover, each state is associated with a field in the 
conformal field theory describing the critical behavior of the 
chain. These conformal fields satisfy certain fusion rules, 
which, generally, are different from the fusion rules of the 
anyons themselves (typically, the number of conformal fields 
differs from the number of types of anyons). As a result, the 
topological symmetry sectors must be associated with the con- 
formal fields in a manner that both the su(2)fc fusion rules of 
the anyons and fusion rules of the conformal fields are satis- 
fied. For the case of the su(2)j. spin- 1/2 anyonic chain, this 
constraint is obeyed for the following reason: the relevant 
critical theories are so-called coset theories, which contain a 
su(2)fc theory and other theories such as u(l). This implies 
that the fields in the critical theory inherit su(2)j. topological 
symmetry labels; thus, the topological symmetry sectors can 
be assigned in a consistent manner. 



A. The ferromagnetic case 

The ferromagnetic su(2)fe spin-1/2 anyon chain is de- 
scribed by the coset theory su(2)k /it (l) 2 /c (details can be 
found in appendix E4i. The fields in this conformal field 



theory carry two labels, a su(2)fe label I, and an u(l) la- 
bel m, where I = 0, 1, . . . , k, m — 0, 1, . . . 2k — 1 and 
I + m = mod 2. Under fusion of two fields with labels 
(li, mi) and (I2, m^), the labels m,\ and 7712 are added mod- 
ulo 2k, while the labels l\ and I2 satisfy the fusion rules of 
su(2)fc. Thus, the fields (l,m) can be assigned a topological 
label I, and this assignment automatically obeys the correct 
fusion rules. 

The momentum quantum numbers of the fields cannot be 
predicted from the conformal field theory itself. Different re- 
alizations of a particular CFT may vary in the assignment of 
momenta to conformal fields. For example, the su(2) 2 spin- 
1/2 chain is described by the Ising CFT for both AFM and FM 
interactions, but the states corresponding to the a, or j = 1 /2 
field occur at different momenta, as illustrated in Figure [TT] 

We first consider the case of k even. The Hilbert space 
is given by labelings of the the fusion chain as displayed in 
Figure [I] where, in the case of the spin-1/2 anyon chain, the 
'incoming' labels are spin-1/2 anyons. As a consequence, 
the labels xi alternate between integer and half-integer values. 
Thus, there are two decoupled 'sectors': In one sector, the la- 
bels of the odd sites correspond to integer-spin anyons, while 
in the other sector, the odd sites correspond to half-integer- 



spin anyons. Because of this, each field in the conformal field 
theory will appear twice in the spectrum, once at momentum 
K and once at momentum K + it. 
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TABLE V: Scaling dimensions in the Z4 parafermion model. 

As was already discussed above, for the case of even fc, the 
topological sector y of a field labeled by (/, m) is determined 
by I, namely y = 1/2. Our numerical results show that the 
momenta of the fields are either given by K(m) — or 
by K(m) = ir + as can be seen in the right hand side 



panel of Figure 18 To establish that the FM spin-1/2 chain 



is stable under perturbations preserve both spatial and topo- 
logical symmetry, we need to show that there are no relevant 
operators with the same momentum and topological quantum 
numbers as the ground state. 

In the case of odd fc, anyon spins j are automorph to anyon 
spins k/2 — j (see Appendix A), and therefore the labels of 
the conformal fields are given by (I, m) where both I and m 
are even. The topological sectors are given by 1/2, and the 
momenta of the fields are given by K(m) = 

From the above discussed relations between field labels 
(I, m) and the quantum numbers (topological sectors and mo- 
menta), it becomes apparent that each momentum and each 
topological sector appears at most once. This implies that the 
critical behavior is indeed stable to perturbations which pre- 
serve both spatial and topological symmetry. 



B. The anti-ferromagnetic case 

In this section, we show that the criticality of the anti- 
ferromagnetic su(2)fc spin-1/2 chain is stable under pertur- 
bations that do not break the symmetries of the model. The 
model is described by the fc-critical Ising model, which 
can be formulated in terms of a coset-model su(2)i x 
su(2)k-i/ su(2)k (some details of this coset model can be 
found in appendix |E l\ . The conformal fields in this CFT are 
labeled by (r, s), where the r label (1 < r < fc) is associated 
with su(2)^_i, while s (1 < s < k + 1) is associated with 
su(2)fc. There is also a label associated with su(2)i, however, 
this label is fixed by the constraint t = r + s mod 2. 

The topological sectors are given by (s — l)/2. Since s is 
the conformal label associated with the denominator su(2)fc of 
the coset, the fusion rules of the coset CFT are consistent with 
the fusion rules associated with the topological sectors. 

In the case of even fc, all fields appear twice in the spectrum, 
(once at momentum K and once at K + ir) as a result of the 
'doubling' of the Hilbert space. Our numerical calculations 
yield the following. The topological sector of each field is 
determined by s, namely y = (s — l)/2. The momentum of a 
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FIG. 19: (color online) The su(2)s spin- 1/2 chain: Energy spectra have been rescaled to match the conformal field theory prediction given 
in Eq. {22). Green squares indicate the location of the primary fields, red circles the descendant fields. The topological symmetry sector is 
indicated by the violet index. Data shown is for system size L = 26 and L — 25, respectively. 
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TABLE VI: Scaling dimensions for the tetra-critical Ising model. 



dimerizing the system, the most relevant field has labels 
(r, s) = (2, 1) and thus a scaling dimension 2h = (k + 
4)/(2(fc + 1)) that is relevant for all k. Therefore, a pertur- 
bation which breaks translational symmetry may open up a 
gap- 



VII. DISCUSSION 



field labeled by (r, s) is given by either K = (r + s mod 2)tt 
or by K = (r + s + 1 mod 2)7r; the system size determines 
which one of the two possibilities occurs ( we verified this 
behavior for k — 2,4). 

For odd k, the association of field labels (r, s) with topo- 
logical and momentum sectors coincides with that for even 
k. However, only odd values of s appear, due to the above 
mentioned automorphism of anyon spins. These results were 
verified for k = 3 in Ref. ([I6i l. and for k = 5 in this study 
(see figure [l~9|>. 

To confirm that the criticality of the AFM spin- 1/2 chains 
is stable under perturbations which preserve the spatial and 
topological symmetries of the model, we have to analyze the 
scaling dimensions of the fields which have the same topolog- 
ical quantum number as the ground state. The ground state 
has label s = 1 (i.e., topological sector y = 0). The scal- 
ing dimensions of the fields with label s = 1 are given by 
2h = (r 2 (fc + 2) - 2(k + l)r + k)/(2{k + 1)), which for 
r > 1 increases monotonically. The most relevant field in the 
same momentum sector as the ground state thus carries the la- 
bels (r, s) = (3, 1), and has scaling dimension 2h = 2 + ^rpj, 
which is irrelevant for k finite, and becomes marginal in the 
limit k — > oo. Again, we conclude that the AFM spin-1/2 
chains are stable with respect to perturbations which preserve 
both topological and translational symmetry. 

When breaking the spatial symmetry of the model by 



The anyonic analogs of the SU(2) Heisenberg spin-1 model 
have a rich structure, as can be seen from the phase diagrams 
of the ordinary bilinear-biquadratic spin-1 model, the generic 
even k > 6 anyonic model, the generic odd k > 5 anyonic 
model, and the special case k — 4 (displayed side by side in 



Fig. [20 1. 

The distinct nature of the phase diagram for k — 4 origi- 
nates in the symmetry of the fusion rules of the su(2)fc theory 
under the exchange j <H- | — j which, for k = 4, maps j = 1 
onto itself. It is also the lowest value of k for which a generic 
fusion rule 1x1 = + 1 + 2 applies (compare with k = 3 
where 1x1 = + 1), thus making it possible to define an 
anyonic spin-1 model. Moreover, the central charge of the 
defining su(2)4 algebra is an integer (c = 2), and the quan- 
tum dimensions of the su(2) 4 anyons are all integers or square 
roots of integers (we discuss the various anyon models in more 
detail in Appendix [A]). We note that fusion models with such 
quantum dimensions typically do not permit 'universal quan- 
tum computation', a pro perty which requires a 'fine tuning' of 
the braid propertie d 45 ! 46 !. Models analogous to the case k = 4 
have been studied from the integrability point of viev J 47 ! 48 ! 

Upon increasing the level k, the su(2)fc anyon model in- 
creasingly resembles the ordinary SU(2) spin algebra. In 
terms of the quantum group language, the limit k — > oo cor- 
responds to q — > 1, where q = e* l ^ k+2 \ For q = 1, the 
quantum group reduces to the ordinary SU(2) algebra. One 
would therefore intuitively expect that the phase diagram of 
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FIG. 20: Phase diagrams of the various spin-1 models considered in the paper: (a) the bilinear-biquadratic spin-1 Heisenberg chain, (b) the 
generic even k > 6 anyonic spin-1 chain, (c) the generic odd k > 5 anyonic spin-1 chain, (d) the special case k — 4. 



the generic k case has the same structure as the phase diagram 
of the SU(2) bilinear-biquadratic spin-1 chain. The numerics 
presented in the paper shows that this is indeed the case for 
both even and odd k, with one notable exception: for even 
k, we find a gapped dimerized phase that is separated from 
the Haldane gapped phase by an extended critical region; in 
contrast, for odd k, we observe only an extended critical re- 
gion but no dimerized phase. The fact that the anyonic spin-1 
models behave differently for even and odd k is very interest- 
ing in its own right. To the best of our knowledge, this is the 
first time that such an 'even-odd' effect in the level k has been 
observed. 

In the following, we discuss some of the differences be- 
tween the cases of odd k versus even k. It is instructive to 
consider the model for the lowest (generic) value of even k, 
i.e., k — 6. From the symmetry of the fusion rules under 
the exchange j «-» | — j it follows that j ' = 1 and j = 2 
are exchanged. This implies a 'symmetry' in the phase di- 
agram of the k = 6 model under exchange of the projec- 
tors and p( 2 \ The parametrization chosen in this paper, 

H = J2i cos — sin^2,i-P ( ^ 1 \ renders the phase dia- 
gram symmetric in the line through the points 6*2,1 = 37r/4 
and 6*2,1 = 77r/4. It is important to realize that this 'symme- 
try' only applies the values of the energies which appear in 
the spectra, but not to the momenta and the degeneracies of 
the energy levels. In particular, the gapped dimerized phase 
is the 'mirror phase' of the gapped Haldane phase. We also 
note that the same mechanism resulted in a symmetric phase 
diagram for k = 4, if plotted in terms of the projectors p(°> 
andP< 2 ). 

For even k = 6, there is an extended critical region between 
these two gapped phases; however, we were not able to deter- 
mine its precise critical behavior. For k — 8, the extend of this 
critical region is smaller, and it is therefore not inconceivable 
that for very large even k, this critical region will shrink to 
a single critical point separating the two gapped phases, as is 
the case for the SU(2) spin-1 bilinear-biquadratic model. 

As indicated by the above terminology, the gapped phase 
around the angle 6*2.1 = is the anyonic analogue of the Hal- 
dane gapped phaseR The ground states at 6*2,1 = can be 
obtained exactly, and th ey are the anyonic analogues of the 
AKLT state 21 . In section HIE 4 we studied this 'AKLT' point 



of our anyonic models with open boundary conditions. We 
obtained edge states similar to the ones observed in the SU(2) 
case. In the case of periodic boundary conditions, we find 
a k + 1-fold degenerate ground state (one for each topolog- 
ical sector), occurring at momentum K = 0. Although the 
Haldane phases of the SU(2) and anyonic models share many 
properties, they differ in their degeneracy for periodic bound- 
ary conditions. Therefore, it is interesting to investigate in 
which way the underlying quantum group symmetry changes 
the classification of gapped phases in one-dimensional spin 
systems^. 

The dimerized gapped phase of the anyon model exhibits 
exactly the same values of the energy levels as the gapped Hal- 
dane phase, as pointed out above. Nevertheless, this gapped 
phase is of a different nature. At 6*2,1 = 37r/2, there is a 
(unique) zero energy state at momentum K = ir, i.e., the 
phase is dimerized - like the corresponding phase in the SU(2) 
spin-1 model. In addition, there is a set of degenerate zero en- 
ergy ground states at momentum K = 0, where the number 
of states depends on k. 

Almost two decades ago, Koo and SaleuPl considered 
a spin-1 type loop model that was based on the 'fused 
Potts model'. The underlying algebra of their model is the 
Birman-Murakami-Wenzl (BMW) algebra, which replaces 
the Temperley-Lieb 49 algebra that appears in the study of the 
Potts model in its loop representation. For details on the 
BMW algebra, we refer to 50 . The Koo-Saleur model contains 
a continuous parameter Q, which is closely related to the dis- 
crete level k in the anyon models we consider (see below). 
More specifically, the model Koo and Saleur consider is 

H KS = ^(Q-l)(sina;-cosw)P J (0) -(g-2)coswP (1) . 

(33) 

The projectors PW and PW project two neighboring spin-1 
loops onto the spin-0 and spin-1 channels, respectively (see^ 
for explicit expressions of these projectors in terms of the 
BMW algebra). The number of Potts states Q is related to 
the quantum dimension of the spin- 1/2 anyons di/2 (or the 
parameter d appearing in the Temperley-Lieb algebra) via 
Q = d 2 1/2 = 4cos(7r/(fc + 2)) 2 , and thus Q = 1,2, 3,4 cor- 
responds to k — l,2,4,oo. In particular, the case Q = 4 
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corresponds to the ordinary SU(2) spin-1 chain. We note that 
the anyonic chains can only be defined for integer k > 4, 
and recall that we parametrized the anyonic spin-1 model as 

H = Ei C0S (°2,i)P- 2) - sin((9 2 ,i)P (1) . By making use of 
the relation 1 = P (0) + + P (2 \ one finds the following 
relation between the parameters of the models 

cos 02,1 = — (Q — l)(sincj — cosw) 

sin 6*2.1 = — cosu + (Q — 1) sinw . (34) 

Despite the similarities between the model of Koo and 
Saleur and our anyonic model, they behave rather differently. 
The phase boundaries between the various phases observed 
in the Koo-Saleur model depend smoothly on the continuous 
parameter Q, while the phase diagrams of the anyonic spin- 
1 models depend on whether k is even or odd. In addition, 
the Koo-Saleur model displays non-unitary critical behavior, 
while the critical behavior of the anyon models is described 
by unitary CFTs. The explanation for this difference in be- 
havior should be sought in the representations used in the two 
models. In the Koo-Saleur model, a representation which es- 
sentially behaves like a SU(2) representation is used (which 
permits to define the model as a function of the continuous 
parameter). In the anyonic version, the truncated su(2)fc rep- 
resentations play a central role. For a related discussion in the 
general context of loop models, we refer td sl l 52 [ 

These observations suggest that a deeper investigation into 
the differences and similarities of the two models is warranted, 
especially because the Koo-Saleur model exhibits various in- 
tegrable points^. One of the integrable points identified in 
l40l corresponds to the supersymmetric critical point forming 
the boundary of the Haldane phase. The location of this in- 
tegrable point, in terms of the parameters used in this paper, 
istan6» 2 ,i = where d x = 1 + 2 cos (2?r/(fc + 2)) 

(see^l). For k > 4, this location depends only weakly on 
k, namely, for k = 4, one obtains #2, 1 = — arctan(3/4) rj 
— 0. 20487T, while in the limit k — > 00, one obtains #2,1 = 
— arctan(2/3) « — 0.18727T. The location of the critical end 
point of the Haldane phase we obtained in this paper are con- 
sistent with the location of this integrable point. 

To solve the anyonic spin-1 chain at this integrable point, 
one approach is to map the model to a fused RSOS model, as 
studied i n 23 * 54 ! (see also 55 56). This subject will be described 
in a separate publication^! 
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Appendix A: su(2)k anyons 

In this appendix, we briefly review of the properties of 
su(2)fc anyons - the building blocks of the anyonic chains con- 
sidered in this paper - for arbitrary level k > 2. We explicitly 
discuss the levels k — 2, ... ,7. For a general discussion of 
anyon models, see e.g. Refs. 1 141 151201 

The anyons of the su(2)fe theories are closely related to 
ordinary SU(2) spin degrees of freedom; thus we label the 
anyons by their 'generalized angular momenta', or simply 
'spin' value j = 0, |, 1, 4 . We note that in the su(2)fc the- 
ory, there is a maximum allowed value of the 'spin', namely 
k/2, a feature not present for ordinary SU(2) spins. 

Ordinary spins can be combined using tensor products. In 
general, combining two spins gives rise to several different 
spins. An analogous phenomenon occurs if we combine two 
anyons of the su(2)fc theory. In the following, we will assume 
that k is fixed, but arbitrary; i.e., the anyons combined belong 
to the same theory. The rules for combining two anyons - also 
denoted as 'fusion rules' - are closely related to the SU(2) 
tensor products, namely 

naa(J+j' ,k—j—j') 
3"=\j—j'\ 

The only difference to the case of ordinary SU(2) spins is the 
cutoff in the upper limit of the sum in Eq. ( |A 1 1 >. The cutoff is 
the result of the finite number of types of anyons in the su(2)fe 
theories. The fusion rules in Eq. ( |A 1 [ > are associative. 

The fusion rules can be represented in terms of the fusion 
matrices Nj which, in the case of su(2)fc anyons, have entries 

(the so-called fusion coefficients) N J - ., = 1 if and only if the 
fusion of labels j and j' gives rise to the label j", and zero 
otherwise. In general, fusion coefficients bigger than one are 
possible, but they do not appear in the context of this paper. 

Fusion is commutative and associative, i.e., fusing several 
anyons in different order gives rise to the same result. This 
implies that the fusion matrices Nj commute and that they 
can be diagonalized simultaneously. Diagonalizing the fusion 
matrices yields the quantum dimensions dj, 

Nj d = djd, (A2) 

where d is a vector whose components are the quantum di- 
mensions dj . The total quantum dimension T> is defined as 




(A3) 



For su(2)fe anyons, the quantum dimensions are given by 

do = 1 , 

di/2 = 2C0S (^) • 

dj = d 1/2 dj_ 1/2 - dj-t, j>l. (A4) 
Explicitly, one obtains 
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where we note that the dimensions dj depend on the level fc, 
which we have suppressed in the notation. 

The matrix which diagonalizes the fusion rules is called the 
modular S-matrix. Its entries for the su(2)^. theories are given 
by 



fc 



(2j + l)(2/ + 1)tt 
k + 2 



(A6) 



For odd k, there exists an automorphism relating anyons 
with spin j to anyons with spin | — j. The automorphism 
thus relates integer and half-integer spins, reducing the study 
of odd-fc anyon systems to only integer (or only half-integer) 
any on spins (it also means that there are only fc/2 distinct 
anyon types for odd k). In this manuscript, we consider 
anyons with integer spin when studying odd-fc systems. 

The Hilbert space of a multi-anyon system is non-local, and 
it can be represented by a a trivalent graph with each line 
segment representing an anyonic degree of freedom. Such a 
graph is called a fusion diagram. The labeling of the segments 
has to be such that the fusion rules are obeyed at all the ver- 
tices. In Figure |2Tj we display the fusion diagram that defines 
the Hilbert space of the models studied in this paper. 

Each distinct labeling of the fusion diagram defines a basis 
state = \xq, X2, Xl-i). We define the basis states \ip) 
to be orthogonal, i.e., the inner product of two basis states is 
one if the labels of the two states are identical, and zero other- 
wise. The number of basis states in a chain of spin-j anyons of 
length L grows asymptotically as d^, where dj is the quantum 
dimension of the anyon of type j. It is important to note that 
dj generally is not an integer, as would be the case for ordinary 
SU(2) spins. This means that it is not possible to associate a 
local Hilbert space with each anyon, and that the total Hilbert 
space is not a simple tensor product of local Hilbert spaces. It 
also implies that there are no 'internal s z quantum numbers' 
in anyonic Hilbert spaces. The reason behind all these fea- 
tures is that the fusion rules enforce non-local constraints on 
the possible labelings of the fusion digrams. 

In order to define Hamiltonians acting on anyonic Hilbert 
spaces, the anyonic analog to the 6j-symbols for ordinary spin 
degrees of freedom has to be considered. The anyonic version 
of the 6j-symbols is the so-called F-transformation, which 
relates the two different ways three anyon spins, ji, j 2 , jz, 
can fuse into a fourth anyon spin j'4. The F-matrix can be 
defined as a result of the associativity of the fusion rules: 



J4. 




(A7) 



In the case of su(2)^,, the F-matrices are uniquely determined 
by a consistency relation, namely the pentagon equation, and 
by imposing unitarity. A useful expression (for general fc) can 
be found in 57 , and is given in the appendix |b| 

A further basis transformation of interest is the so-called 
S-transformation which relates the 'flux' of anyon spin j 
through a loop of anyon spin I to the case without anyon loop 



3 3 3 3 3 



X Xi X 2 X3 ... 

FIG. 21: Basis (fusion diagram) of a chain of spin-j anyons (j = 1/2 
in the case of the su(2)^ spin-1/2 chain, discussed in section |vT|and 
j = 1 in the case of the su(2)k spin-1 chain, discussed in section|Hr}, 



by 



&i,0 



(A8) 



The matrix elements of this transformation are the elements 
of the modular S'-matrix 14 1S . 

In the following, we give matrix representations of some of 
the above discussed properties of a model of su(2)^ anyons. 
Upper indices in round brackets denote the level fc. 
Level k = 2: Ising anyons - This class of anyons comprises 
the spin-0 anyon, the Ising anyon (spin-1/2) with non-Abelian 
braiding properties, and the fermion (spin-1). The non-trivial 
fusion rules are given by 



1x1=0+1 
2 2 



1 



x 1 



1x1 = 0. 



The corresponding quantum dimensions are given by 



7(2) 



,(2) 



1 



d[ 2) 



(A9) 



(A10) 



and the S'-matrix takes the explicit form (the entries are or- 
dered according to ascending value of the anyon spins) 




(All) 



Level k = 3: Fibonacci anyons - This class of non-Abelian 
anyons exhibits only two distinct particles, with spins and 
1 (the Fibonacci anyon) respectively (the spins | and | are 
automorph to spins 1 and 0, respectively). Thus, there is only 
one non-trivial fusion rule, 



1x1=0+1 
and the quantum dimensions are given by 



(A12) 



7 (3) 



1 



<4 3) - (1 + \/5)/2 



(A13) 



Using the notation </> = (! + v5)/2, the S-matrix reads 



1 



1 
\/2T0 _1 



(A14) 
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Level k = 4 — The k = 4 anyon model contains five anyon The quantum dimensions can be obtained from Eq. ( A4 1, 
spins, namely j = 0, |, 1, §, 2. The fusion rules are given by 



1 



+ 1 



i + « 

2^2 

0+1 + 2 



3 ? 

i 

1 2 
(J 



(A15) 



The (non-trivial) quantum dimensions can be obtained from 
Eq. <|A4]i, 



4 4) = 4 4) = i 



2 2 

Finally, the 5-matrix takes the form 



5^ = 



(A16) 



2V^5 
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V3 




V i 


-V3 
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-73 


1 / 



(A17) 



Level k = 5 - This class of non-Abelian anyons gives rise to 
three distinct anyon particles with spins 0, 1 and 2 (which are 
automorph to the spins |, I and i, respectively). The non- 
trivial fusion rules are given by 

1x1 = + 1 + 2, 1x2 = 1 + 2, 2x2 = + 1. 

(A18) 



The quantum dimensions take the following values, 

d (5) - 1 

4 5) = (4 5) ) 2 - 1 = 1 + 2cos(2tt/7) , 
4 5) = 2cos(7r/7) . 

The S'-matrix of the su(2)g theory is given by 



S ■■ 



D(5) 



1 

(5) 
,(5) 



4 5) 


4 5) 


-4 5) 




i 


-4 5) 



(A19) 



(A20) 



where denotes the total quantum dimension — 

\Jl + {df ] ) 2 + (4 5) ) 2 of the su(2) 5 theory (restricted to the 
integer 'spins'). 

Level k = 6 - The anyon model with k = 6 has seven anyons 
labeled by j = 0, h, 1, 1, 2, |, 3. The fusion rules read 
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4 6) = 4 6) = i 



2 2 

4 6 ' = 4 6) = 1 + ^2 4 6) = vW 2 + V2. (A22) 



The entries of the 5-matrix are given by Sij 

ySsin( ^ 2i +gg± 1 > ), for i, j = 0, 1/2, 1, . . . , fc/2, 
fc = 6. 



Leve/ k = 7 - Finally, we provide some details of the k = 7 
model, which contains four distinct anyons with spins 0, 1, 2 
and 3. The fusion rules are 



x 
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2 3 
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-2 1+2+3 2+ 
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0+1+2+3 1+ 


2 


3 






+ 
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and the quantum dimensions are given by 



4 7) = 1 

7(7) 



4 7) = 1 + 2cos(2tt/9) 



4' ; = 1 + 2cos(tt/9) , 4 ?) = 2cos(tt/9) . (A24) 



The entries of the 5-matrix are given by Sij 
fi»; ,'_n 1 (h_ 1\ / 

(fc+2) 



fc +2 sin( (2z + 1 ^+ 1) " ), for i, j =0,1,..., (k- l)/2, with 



= 7. 



Appendix B: F-matrices of the su(2) fc theories 

In this section, we give an explicit expression for the 
F-symbols, following Ref. 57 We begin with some prelim- 
inary notation. The e-numbers are defined as [nj = 

^^gT -1 — 9 l~ g _i ■ The q-factorials are defined 

q2~q 2 

as [n\\ — [n\ [n — lj • • • |_ 1 J , for integer n > 0, and 
[0J ! = 1. The labels of the anyons a, b, . . . take the values 
0, 1/2, 1, The quantum dimensions are dj = [2j + lj = 



\ fe+2 



A(a,b, c) 



I sm (^feT2 J = dk/2-j- Moreover, we define 



[a + b - cj ! [a - b + cj ! [-a + b + cj ! 



[a + 6 + c + lj ! 



(Bl) 

(A21) where a < 6+c, & < a+c, c < a+6 and a+6+c = mod 1. 
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Using the above introduced notation, the ^-symbols can be written as followP 

(FS% = (-l) a+b+c+d A(a, b, e)A(c, d, e)A(b, c, f)A(a, d, /)vW+TJVW+l] 



E' 



(-l) n [n + lj! 



[a + b + c+ d-n\l[a + c + e + f -n\\\b + d + e + f - raj ! 

1 

d-/J! ' 



(B2) 



[_ra — a — 6 — ej ! [_ra — c — d — ej ! [ra — 6 — c — /J!|_ra — a 
where the sum over n runs over (non-negative) integers such that 

max(a + b + e,c + d + e,b + c + f, a + d + f) < n < min(a + b + c + d, a + c + e + /, b + d + e + f) 
which guarantees that the arguments of the g-factorials are non-ne gative integers. 



Efe 1 



J:) 




FIG. 22: Basis transformation used to obtain the fusion product of 
two neighboring spin-j anyons. 



Appendix C: Microscopic models 
1. Basis and Hamiltonian 

We consider a chain of spin-j anyons, using the basis dis- 
played in Fig. [21] We fix the spin-j anyon to be either 
a spin- 1/2 anyon, or a spin-1 anyon; however, the Hamil- 
tonian defined below can be generalized to any value j £ 
{0, 1/2, fc/2}. Throughout most of this paper, we apply 
periodic boundary conditions, i. <?., xl = xo, where L de- 
notes the number of anyonic quasiparticles in the chain. 

We consider interactions between nearest neighboring spin- 
j anyons. In the case of j = 1/2 (the 'su(2)fc spin-1/2 chain'), 
two neighboring spin- 1/2 anyons may fuse into a spin-0 or a 
spin-1 anyon. In contrast, for the case of j = 1 (the 'su(2)j, 
spin-1 chain'), two neighboring spin-1 anyons may fuse into 
a spin-0, a spin-1, or a spin-2 anyon (for k > 4). In order 
to obtain the fusion product of two nearest-neighbor spin-j 



anyons in the basis shown in Fig. 21 an ^-transformation has 



to be performed, as illustrated in Fig. 22 Consequently, the 
projector onto a particular fusion channel I is composed of two 
F-transformations. This projector, denoted by P i penalizes 



the fusion of anyons at positions i and i 
and it is defined as follows, 



P» |*^0) ••*) 1 ? 2*£ i ^z+1 ? • ■ • i %L — 



E 



1 into an I -anyon, 

_i,ar-,ar i+ i, ...,x L - X ) 
(CI) 



We note that this definition utilizes that F = F^ 1 for su(2)fc. 
The explicit form of the local projectors for the systems stud- 
ied in this paper is given in appendix IS] The Hamiltonians 



discussed in the following section are composed of the sum of 
the local projection operators onto the fusion product I of 
two nearest neighbor spin-j anyons 



2. Hamiltonian of the su(2)fe spin-1/2 chain 



The Hamiltonian of the su(2)fe spin-1/2 chain is given by 



h = jJ2p?\ 



(C2) 



where the projector P^ is defined in Eq. dClll (note that 



I = and j = |), and the coupling constant takes the values 
J = ±1. In analogy to the 'ordinary' Heisenberg spin-1/2 
chain, we denote the case J = —1 in Hamiltonian ( |C2| i as 
antiferromagnetic (AFM) coupling while J = 1 is ferromag- 
netic (FM) coupling. 

In the following, we present matrix representations of the 
Hamiltonians of the spin-1/2 anyon chains for k = 2,3, 4, 5. 
The matrix formulation for the su(2)3 spin-1/2 chain was 
first introduced in 16 (see also 58 ). Local basis elements are 
labeled by a;,, where Xi £ {0, 1/2, 1, fc/2}. The or- 
der of anyon spins in the matrix representation is ascending. 
We also introduce the operators n\ acting on local state \xi): 
n\\xi) — e\xi) where the eigenvalue e = 1 if the local basis 
element Xi = j, and e = otherwise. 

From the definition of the summands of the Hamiltonian 
(eq |Cl| l, it is apparent that the matrix representation of the 
projector depends on the basis elements Xi-\, xi and Xj+i. 
In particular, non-trivial contributions to Pf exist only for cer- 
tain values of Xi-i and Xi+i, and thus each contribution to the 
projector will be proportional to nl^n? , v for some values of 
j and f . By specifying both j and f, the possible values of X4 
are fixed by the fusion rules. If there is only one value Xi can 
take (for given j and j'), we omit the identity operator that is 
applied to basis element |xj). If there is more than one pos- 
sible value of Xi, we specify the matrix assigning the correct 
energies. 

In the case of even-fc spin-1/2 chains, the fusion rules 
Eq. (A9l impose that the values of the local basis elements 
Xi alternate between integer and half-integer values. For the 
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odd-k spin chains (both spin-1/2 and spin-1 chains), we only even-fc spin-1 chains, the Hilbert space splits into two disjoint 
consider the integer anyon spin subspace (recall the automor- sectors, the integer sector (all Xi take integer values), and the 
phism that applies to odd k anyons, see appendix [A}. For the half-integer sector (all Xj, assume half-integer values). 



a. su( 2)2 spin-1/2 chain 
The Hamiltonian eq. ( |C2[ > takes a rather simple form in the case of su(2)2, namely 

ff (*= 2 ) = J^nUnU + n\_ x nU + +\n\L\n\Z (\ ~A . (C3) 



b. su(2)i spin-1/2 chain 
The Hamiltonian for the fc = 3 spin-1/2 chain is given by 



-Vd 



(C4) 



i 

where d=di = (l + \/5)/2. 

c. su(2) 4 spin-1/2 chain 

In the case fc = 4, the local basis elements alternate between integer spin, x$ € {0, 1, 2} and half integer spin, x.- L+ i g 
{1/2, 3/2}. The Hamiltonian takes the following form 

rx(fc=4) _ jr„0 _0 ,_2 „2 , 1 1 1 fl -l\ 1 1/2 1/2 / 1 -\/2\ 1 3/2 3/2 / 2 

w — J ,7^ + 2 1+1 \-l 1 / 3 +1 1-^2 2 / 3 i_1 1+1 \-\/2 1 / 

(C5) 



d. su( 2 )s spin- 1/2 chain 
Using the notation d\ = 1 + 2 cos(27r/7) and d 2 = 2 cos(7r/7), the Hamiltonian reads 



H<-> . J^nUnl, + ^nUnU (_ ^ + "f^ (CO 

3. Hamiltonian of the su(2) fc spin-1 chain 

We define the Hamiltonian of the su(2)fc spin-1 chain as follows, 



The projectors P^ and P^ are defined in Eq. (CI I, where I = 1 and I = 2, respectively. This Hamiltonian is the su(2)fc 



anyonic equivalent of the bilinear-biquadratic spin-1 chain. Throughout the paper, we parametrize the Hamiltonian by the angle 
9 as follows: J\ = — sin.(#2,i)> J2 — cos(#2,i). 

The Hamiltonian Eq. ( |C7[ ) is defined for levels k > 4, in which case the fusion of two spin-1 anyons may result in a spin-0, a 
spin-1 or a spin-2 anyon (for level k — 3, spins | and 1 are automorph, i. e., the spin-1 chain is equivalent to the spin-1/2 chain; 
moreover, the fusion rules imply 1x1 = + 1. For k = 2, the fusion of two spin-1 particles is trivial, 1x1 = 0). 
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a. The su(2)i spin-1 chain 



We shall now present a matrix representation of the Hamiltonian of the su(2)4 spin-1 chain, using the same notation as in the 
previous subsection. In the integer sector (IS), the projectors onto the different channels can be written as follows, 



Pi% = «?-!»?+! + + \n\-in\ +1 ^-^2 




1 



i o -r 



p lis = n i-i n i+i + n l-i n °i+i + + n t-i n l+i + ~ n \-i n \+i 

'10 1 



(2) 
,IS 



1 y/2 1 



1 



1 y/2 1 

In the half-integer sector (HIS), we can write the projectors as follows, 



P (0) --(n 1/2 n 1/2 +n /2 n /2 )( 1 ^ 



>« - l^V2„l/2 , 3/2 3/2 v A 1 



i,fflS = + (I [) + + n i->i+l) 

d(2) 1 / 1/2 3/2 . 3/2 1/2, A 1 

^,HIS = o( n <-l n <+l + n i-l"i+l) I 1 1 



3/2 1/2 x 



1 -1 
-1 1 



(C8) 



(C9) 



(CIO) 



(Cll) 
(C12) 
(C13) 



b. The su(2) 5 spin-1 chain 
Using notation d\ = 1 + 2 cos(27r/7) and d 2 = 2 cos(7r/7), the projectors are given by 
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(C14) 



,(2) 



/ d 2 d 2 d 3/2 d 2 di( j 
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(C15) 



c. The su(2)§ spin-1 chain 

In the following, we use the notation di/ 2 = 2cos(7r/8), d\ = 1 + 2cos(7r/4) = 1 + \/2 and d 3 / 2 = 2\/2cos(7r/8). In the 
integer sector (IS), the projectors onto the different channels can be written as follows, 



P tJS = n t-l n i+l + n l-l n "+l + n i-l n t+l + n l-l n i+l + ^K'-l^i+l + 



1 -1 
-1 1 
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n i-i n i+i 
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+ 



A A A ~ U i-l n i+l 
«l/2«l«3/2 

°-i n i+i + n \-i n \+i + n i-i n i+i + n i-i n l+i + \( n l-i n t+i + n 1-i n \+i) (l \ 

^1 /2^3/2 
1^1/2^3/ 
\fd^l/2<h/2 



p( 2 ) - „ 



d\ d 1 -d ly /d 1/2 d 3/2 \ 

di J -y/d 1/2 d 3/2 (C16) 

di^/di/ 2 d 3 / 2 -^d x / 2 d 3 / 2 di/ 2 d 3 / 2 J . 



d\/ 2 did 3 / 2 
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In the half-integer sector (HIS), the projectors onto the different channels can be written as follows, 
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y / di/ 2 did 3 / 2 
d. The su(2)r spin-1 chain 

In the following, we use the notation rfi = 1 + 2 cos(27r/9), d 2 = 1 + 2 cos(7r/9) and d 3 — 2 cos(7r/9). The projector 
takes the form 
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The projector P- is given by 
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Appendix D: Exact form of the AKLT states 



For I = p, this results in 



In this section, we present the explicit form of the zero en- 
ergy ground states of the periodic anyonic spin- 1 chains for ft 
odd at the anyonic equivalent of the AKLT point. In the main 
text, we discussed the case k = 5. 

At the AKLT point, the Hamiltonian contains only the pro- 
jector onto the anyon spin-2 channel, i.e., the fusion of neigh- 
boring spin-1 anyons into a spin-2 anyon is penalized. First, 
we note that the fusion of anyons of types 1 and (ft — l)/2 (the 
latter being the largest integer'spin' for an anyon in the su(2)fc 
theory) results in 1 x (ft - l)/2 = (ft - 3)/2 + (ft - l)/2. 
In addition, we have that 1 x (ft - 3)/2 = (k — 5)/2 + 
(ft — 3)/2 + (ft — l)/2. Thus, a local basis for which 
Xi-i = Xi+i = (ft — l)/2 implies that Xi — (ft — 3)/2 or 
Xi = (ft— 1)/2. It follows that after the local basis transforma- 
tion, Xi can only take two possible values, namely Xi = 0, 1 
(consider 2 x (ft - l)/2 = (ft - 5)/2 + (ft - 3)/2, and let 
= (ft — l)/2, then x i+ i could only take values 
(ft - 5)/2 and (ft - 3)/2 but not (ft - l)/2). This, in turn, 
means that a choice of local variables Xi-\ and Xi+i does not 
give rise to non-zero contributions at the AKLT point as fu- 
sion of neighboring spin-1 anyons in the chain cannot result 
in ii — 2. We thus obtain a zero-energy ground state of the 
form \vo) = \(k- l)/2, (ft - l)/2, . . . , (ft - l)/2). 

To construct the other ground states, we make use of the 
topological symmetry operators Y\ . These operators mutually 
commute, and they commute with the Hamiltonian. The state 
vq is not an eigenstate of the operators Y\ (with I > 0), and 
hence alternative zero-energy ground states of the Hamilto- 
nian are given by \vi) = Y\ \ vq). These ground states \v{) 
(/ > 0) can be obtained explicitly. The local basis states take 
values Xi = p— I or Xi = p—l + 1, where p = (ft — l)/2. The 
states \vi) are a sum over all possible labelings of the fusion 
tree with these two values of X4. We introduce the follow- 
ing notation: #1 denotes the number of local basis states for 
whicha;i = I, and m) denotes the number of local basis 
states for which Xi — I and Xi + i = m, where we use periodic 
boundary conditions, xl = xq. For I > 0, we obtain 



Xi£{p-l,p-l + l} 



fl({Xi}) \x ,xi, 



,X L -i) 



(Dl) 



The coefficients fi{{xi}) (0 < I < p) are given by 



d„ 



#(p-l+i, P -i) 



J dp-idp-i+i p 
di+i. 



fp({Xi}) 



ip-i 



L/2 



dp-i 

d p y/di' 



\#1 



We denoted the ground states by \vi) for the following rea- 



son. In section II A we explained that the operators Yi can be 
thought of as fusing an anyon with 'spin' I into the chain, ef- 
fectively changing the 'overall fusion channel' , or flux thought 
the chain. If we take a state \iij 2 ), and act on it with the oper- 



ator Yj lt we find that Yj t \vj 2 ) 



E 



i ) , where the 



sum is over those j'3 which appear in the fusion ji x J2- Thus, 
the ground states of the AKLT anyonic spin chain form a 'rep- 
resentation' of the fusion algebra of su(2)fc. This implies that 
eigenstates of the topological operators Y[ can be constructed 
because the modular S-matrix diagonalizes the fusion rules. 
In particular orthogonal (not normalized) ground states at the 

AKLT point are written as IV'AKmy) = Y^^=q^ 2 

where Sij is the modular S matrix for su(2)fc, and the sum is 

over integer values. 



Appendix E: Conformal field theories of interest 

In this appendix, we summarize the most important aspects 
of the conformal field theories relevant to this paper. In the 
following, 'primary fields' refers to Virasoro primary fields. 
Detailed discussions of conformal field theories can be found 
in Ref. |59]and Ref|60] 



1. Virasoro minimal models 



The unitary minimal models^, which can also be described 
in terms of the coset ^"( 2 )i X;i "( 2 )fc-i nave a cen ]- ra i charge 

su(2) k ' & 

c = 1 — re , Tvrr, i n \ (ft > 2). The primary fields are labeled 



(fe+l)(fe+2) 

by integers r and s, where 1 < r < ft and 1 < s < ft + 1. 
Their conformal dimensions are given by 



(r(fc + 2) - (ft + l).s) 2 - 1 
4(ft + l)(ft + 2) 



(El) 



#0-/+x,p-i+i) 



Typically, the minimal models are labeled by a parameter m = 
k + 1. 

Apart from the so-called diagonal models, there exist mod- 
ular invariants that give rise to conformal field theories with a 
different field content 34 35 . More information on these modu- 
lar invariants can be found in table IVIll 



35 



P. p' 


(A,D) 


(r,s) 


mult. 2 


p' = 2(2n + 1) 


(D p '/2+i, Ap_i) 


r odd 


r = p' — r 


p = 2(2n + l) 


(y4 p /_ 1 , £) p /2+l) 


s odd 


s = p — s 



TABLE VII: Modular invariants of the Virasoro minimal models. For 
a given pair of indices, (p,p') = (m + 1, m) and n integer, only 
the fields with indices (r, s) as specified in the third column appear 
(1 < r < p', 1 < s < p). Some fields have multiplicity two, as 
indicated in column four. 



2. Af — 1 superconformal minimal models 

The Af = 1 superconformal minimal models^ are de- 
scribed by the coset 

su(2) 2 x su(2) fe _2 



su(2) k 



and have central charge c 
have conformal dimension 

(r(k + 2)-sk) 2 



12 



k(k+2) 



The primary fields 



1 



(1 - (-l) r - s ) , (E2) 
1. The fields with 



{r ' s> 8fc(fc + 2) ' 32 1 

where 1 < r < k — 1 and 1 < s < k 
r + s even, i.e. the fields in the Neveu-Schwarz sector, have a 
super partner, whose conformal dimensions are given by 



for r + s even . (E3) 



4. The Z fc parafermion CFT. 

The Zk parafermionJSl can be described in terms of the 
coset 

su(2) k 
u(l)2fe ' 

where u(l)2fe denotes the c = 1 boson, compactified on a 
circle of radius R = v2fc. The central charge is given by 



2(fc-l) 



, and the conformal dimensions of the primary 



fields are given by 



l (l,m) 



1(1 + 2) 

4(fc + 2) 



m 

~4k 



(El) 



Here, the indices run over values I 
-1 + 2,-1 + 4,. ..,1. 



0,1,. 



, k and m 



5. The Z2 orbifold theories 

We will briefly discuss the Z 2 orbifold of the compactified 
boson at squared radius R 2 = 2p. For a detailed account, 
we refer td=4 The number of primary fields is given p + 7, 
where p — 1,2,.... For p = 1, the CFT is Abelian, and 
it is equivalent to a the compactified boson theory with eight 
primary fields. In general, the following fields are present. 



3. S3 minimal models 



The class of S3 symmetric minimal are de- 

scribed by the coset theory 

su(2)4 x su(2)k-4 



su{2) k 



and have central charge 

c = 2 



24 



(k-2)(k 



(E4) 



12 



There are two main sets of primary fields. The first set has 
conformal weights 

_ (r(k + 2)- s(k - 2)) 2 - 16 1 - cos 4 (7r(r - s) /4) 
{r ' s) ~ 16(fc-2)(fc + 2) + 

The second set has scaling dimensions 

h '(r,s) = h (r,s) + 

1 + sin 2 (7r(r - s)/4) 



(E5) 



+ $r,i5 s ,i + 8 r i5 s o + 2S r . 2 5 s 



(E6) 



where for both sets 1 < r < k — 3 and 1 < s < k + 1. 
There are additional (Virasoro) primary fields, with scaling 
dimensions differing by integers from the scaling dimensions 
listed above. These additional primary fields are not relevant 
to this work. 



• The identity field 1, with scaling dimension hi = and 
quantum dimension d\ = 1. 

• The field 9, with dimension = 1 and quantum di- 
mension d@ = 1. 

• Two 'degenerate' fields $ x and $ 2 , with scaling dimen- 
sion h$, = I and quantum dimension d$ = 1. 

• The twist fields a 1 , a 2 and r , r 2 , with scaling dimen- 
sions hp = jg and h T — and quantum dimensions 
d a = d T = y/p. 



,p — 1, with scaling 
dimensions h\ = and quantum dimensions d\ = 2. 



• The fields <j)\, with A = 1,2, 

' A 



a. The S-matrix 

To verify that the assignment of the topological symmetry 
sectors of states of the critical su(2) 4 spin-1 anyonic chains 
are compatible with the fusion rules of the orbifold CFTs de- 
scribing the critical behavior, we need the fusion rules of the 
orbifold CFTs. We will not give these fusions rules explic- 
itly here, but specify the modular S-matrix. The fusion rules 
can be obtained from the modular S-matrix by means of the 
Verlinde formula^. 

The modular S-matrix can be written in a compact way as 
follows. In the basis (1, 9, er\ t, <f)\) for the rows and 



36 



(1, 9, cr J , r J , M ) for the columns, where i, j = 1,2 and 
X, fJ, = 1,2, . . . ,p — 1, the modular S'-matrix is given by 



S = 



1 



/ii i Vp Vp 2 \ 

11 1-^-^2 
1 1 (-I)p 6 4J b itj (-1)"2 



Vp ~Vp Ki 



J i,3 



''■J 



\ 2 2 (-1) A 2 
Here, the matrices a, 6, c have the elements 

= v/p72(l + (2^ - l))e-^ 2 








(E8) 



h, „ = (-lf+^ype^/ 2 , 
4cos(7rA/z/p) . 



We note that we used a simplified notation in the above defi- 
nition: for the matrix elements that do not depend i or j, the 
particular element does not depend on i and j. 
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